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Abstract 

In the asymptotic safety programme for quantum gravity, it is important to go beyond 
polynomial truncations. Three such approximations have been derived where the restriction is 
only to a general function f(R) of the curvature R > 0. We confront these with the requirement 
that a fixed point solution be smooth and exist for all R > 0. Singularities induced by cutoff 
choices force the earlier versions to have no such solutions. However, we show that the most 
recent version has a number of lines of fixed points, each supporting a continuous spectrum 
of eigen-perturbations. We uncover and analyse the first five such lines. Sensible fixed point 
behaviour may be achieved if one consistently incorporates geometry/topology change. As an 
exploratory example, we analyse the equations analytically continued to R < 0, however we 
now find only partial solutions. We show how these results are always consistent with, and to 
some extent can be predicted from, a straightforward analysis of the constraints inherent in the 
equations. 
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1 Introduction 



One attempted route to a quantum theory of gravity is through the asymptotic safety programme 
[TJ|6]. Although quantum gravity based on the Einstein-Hilbert action is plagued by ultraviolet 
infinities that are perturbatively non-renormalisable (implying the need for an infinite number of 
coupling constants), a sensible theory of quantum gravity might be recovered if there exists a 
suitable ultraviolet fixed point [l]. 

The task is not just that of searching for an ultraviolet fixed point. 'Conformal gravity', based 
on the square of the Weyl tensor, is perturbatively renormalisable, and corresponds to a Gaussian 
ultraviolet fixed point Vn, for which the Einstein-Hilbert term is a relevant direction. It is apparently 
not suitable however, because the theory is not unitary. 

Suitable fixed points, if they exist, have to be non-perturbative, and unitarity is surely one of 
the more subtle properties that would need to be verified [8j[l6] . There are also phenomenological 
constraints, for example it has to allow a renormalised trajectory with classical-like behaviour in 
the infrared, as emphasised in (9|[l0] . Of particular relevance to the current paper is that there 
should be a fixed point with a finite number of relevant directions (otherwise it would be no more 
predictive than the perturbatively defined theory). Preferably the theory should have only one 
fixed point, at least only a finite number (otherwise again we lose predictivity) . 

Of course the prize of finally combining quantum mechanics and general relativity into a com- 
plete predictive and unified theory cannot be overemphasised, and justify the considerable research 
already made into the possibility of asymptotic safety |3j|6j. It is also the case that functional 
renormalisation groupQ methods provide the ideal framework to investigate this possibility and by 
now provide powerful tools to attack in a different way other longstanding problems of quantum 
gravity that have proved so mysterious and intractable in the past (for progress see [3j|6]). 

It is not possible to solve the full renormalisation group equations exactly however. In a situation 
such as this, where there are no useful small parameters, one can only proceed by considering model 
approximations which truncate drastically the infinite dimensional theory space. At least so far in 
such approximations, scheme independence (i.e. independence on choice of cutoff) and the modified 
BRS Ward identities which encode diffeomorphism invariance for the quantum field [2j , cannot then 
be recovered. Furthermore, aside from the three exceptions which are the subject of this paper, and 



a very recent conformal truncation study in three dimensions 12 , only a finite number of operators 
are kept. Typically these are polynomial truncations, i.e. where everything is discarded except 



1 From here on we will use Wilson's own terminology and call this the "exact renormalisation group" jllj . 
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powers of some suitable local operators up to some maximum degree. Nevertheless, confidence in 
the asymptotic safety scenario has grown strongly: one finds from these approximate models that 
qualitative features persist, and numerical values for universal quantities (renormalisation group 
eigenvalues) are reasonably stable, across many different choices of sets of operators, cutoffs and 
gauge fixings [2]-[6]. 

This sounds very encouraging, however past studies in much simpler field theories have shown 
that it is easy to be misled, by even very high order polynomial truncations, into thinking that fixed 
points exist |13| which can then be shown to be absent in a fuller description of the theory [19] . 

Another example may help to underline this point: in ref. [14] , one of us analysed the exact 
renormalisation group for three-dimensional U(l) theory in the general f(F^ u ) approximation. 
Including such higher powers of F^ u makes the theory perturbatively non-renormalisable. However, 
if any such fixed points existed, one would be able to construct new non-perturbative field theories 
about them. Similarities with the present case should not be lost on the reader. Despite the fact 
that such non-trivial fixed points do exist in polynomial truncations of these equations 14 , it was 
proven that no sensible non-trivial fixed point solution for f(F^ u ) existsj^] 

It is therefore necessary to go beyond polynomial truncations of the exact renormalisation group, 
in order to have confidence that the asymptotic safety scenario is realised. A natural starting point 
in the present case is to keep all powers of the scalar curvature and thus truncate the Wilsonian 



effective action to 15-17 



d 4 x^gf(R). (1.1) 

We will refer to this as the u f(R) approximation" . As Benedetti and Caravelli [17] have emphasised, 
it is as close as one can get to the Local Potential Approximation (LPA) |18[|19| , a successful 
approximation for scalar field theory in which only a general potential V(<£>) is kept for a scalar 
field if (see e.g. pip]). 



The LPA can be viewed as the start of a systematic derivative expansion 20 , in which case this 
lowest order corresponds to regarding the field ip as constant. In rough analogy, an approximation 



of form (1.1) may be derived by working on a four-sphere where the curvature R > is constant, 

although this also has the effect of collapsing more general dependence on R^uXa and R^ to 

this form (since on a maximally symmetric space these are equal to R times products of the 

metric). Using somewhat different methods, three versions of this type of f(R) approximation 

2 There is only the trivial Gaussian solution where / is linear. As we review below, any sensible solution must be 
smooth and defined over the full domain of its argument, here F^ v > 0. 
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have been derived 15-17 . In all these cases the result is a non- linear partial differential equation 
which governs the evolution of f(R,t) with respect to renormalisation group 'time' t = ln(fe//i). 
(Here k is the cutoff value and [i is the arbitrary physical renormalisation scale.) For fixed points 
f(R,t) = f(R), they each reduce to third order non-linear ordinary differential equations. 

This paper is devoted to analysing these equations and perturbations around their solutions, 
adapting techniques that have proved successful in scalar field theory |19| - [23| . Despite the very 
different context, we will gain substantial insight by using these techniques. Indeed, our main 
conclusion is that our simple parameter counting argument remains a powerful way to determine 
what to expect from these equations in terms of fixed points and perturbations, and thus provide 
a first characterisation of the field theories that can be developed around them. This is particu- 
larly important given that in ref. [17], and very recently also in the three dimensional conformal 



truncation study of ref. 12 , partial solutions are apparently found that disobey these simple rules, 
in particular resulting apparently in partial solutions where none are expected and a continuum of 
partial solutions where only a discrete set are expected. Concentrating on the equations of ref. |17| 
we will show that after sufficient care, often requiring very high numerical precision, the parameter 
counting arguments actually do always work as expected. 

Making the appropriate adaptations, the basic argument is as follows (the details are laid out in 
the following sections). In the neighbourhood of a generic point R, the fixed point equation for f(R), 
being a third order ordinary differential equation, has a continuous three-dimensional parameter- 
space of solutions. However as a consequence of the particular non-linearity of these equations, 
almost all of these solutions end at singularities. True fixed points can only be uncovered by insisting 
on the physical requirement that these solutions exist and are smooth^] (i.e. infinitely differentiable) 
for all R > 0. In this case, any so-called fixed singular point of the ordinary differential equation 
imposes a constraint that reduces the parameter space dimension. Furthermore the asymptotic 
behaviour of the solution f(R) for large R can be found analytically and this provides further 
constraints. 

Unless a miracle occurs, which would imply some hidden regularity or symmetry in the equa- 
tions, all these constraints act independently. It is then a simple matter to count them up and 
figure out whether we should expect any smooth solutions to exist for all R, and if so, whether 
they are discrete or continuous in number. Once we know what to expect we can also use these 



In fact to accommodate Litim's cutoff 



24 



, we will need to relax this at certain critical points. This does not alter 
the analysis of whether fixed points exist or not, so we skip over this for now and incorporate it in the full analysis 
later. 
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insights to underpin a careful numerical analysis of these equations and as we will see, confirm 
these expectations. 

What if quantum gravity exactly treated, does indeed have powerful hidden regularities or 
symmetries? We can be confident that if so, they would not survive the approximations we described 
above or indeed others specific to the derivations of the three f(R) approximations 15-17] (again 
see later). 

We will see that we encounter only the generic case for the singular points where a single pole 
reduces the parameter space dimension by one. We will see that the two earliest versions of the 



f(R) approximation [15,16 have four singular points. The equations are over constrained. Three 
singular points are already enough to eliminate all the parameters and leave us with at most a 
discrete set of fixed point solutions. It is therefore not credible that these equations have global 
solutions. 

The origin of some of these singular points are single poles in R due to unphysical divergences 
that occur for the cutoffs at these points (as has been noted in [17]). If ultraviolet fixed points 
do exist in this truncation, either they cannot be reached with these choices of cutoff, or such 
fixed points are formally described by certain singular solutions f(R). Either way, we cannot make 
further progress mathematically. 

A significant advance was made in ref. [17] by choosing the cutoffs so as to avoid these singular- 
ities, however at the expense of not integrating out the lowest scalar mode on the sphere. Together 
with the expected singular point at the origin R = 0, this results in one more singular point in the 
applicable range, at R = R + > 0. 

In the LPA approximation for scalar field theory, an asymptotic expansion can be found for 
V(ip) at large tp, and one finds that number of parameters is again less than would naively be 



expected [19,20 , although this can be understood by considering small perturbations around it. 
The situation here is more involved but we show how the full asymptotic expansion can again be 
understood in this way and we show how to derive it by breaking the equation down into smaller 
parts according to the order of asymptotic behaviour. However here we encounter a surprise. 



Contrary to the conclusions of ref 17 , we find that the full asymptotic series solution for f(R) 
at large R, has three parameters, thus providing no restriction on the dimension of the parameter 
space. Instead these parameters are restricted only to a three-dimensional domain (the inside of a 
cone) . 

Since only the two singular points reduce the dimension of the parameter space, we expect 
one-parameter sets, i.e. lines, of smooth fixed point solutions defined for all R > 0. After some 
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extensive numerical searching, this is exactly what we find. In fact we find five separate lines of 
fixed point solutions, and have evidence for at least one more. 

So far this would be interesting if somewhat problematic for the asymptotic safety scenario. 
However, the parameter-counting is much the same for the eigen-perturbations and it shows that 
the result is unacceptable. The eigen-perturbations obey a linear third-order differential equation 
for any given renormalisation group eigenvalue A. The singular points at R = and R = R + 
each continue to reduce the parameter space by one, whilst the asymptotic series solution does not 
reduce the dimension. Linearity allows a further condition (normalisation) which thus reduces the 
dimension of the parameter space to zero. (If we include A in the parameter-count, we are left with 
one dimension.) In the LPA approximation for scalar field theory, the extra property of linearity 
over constrains the equations thus ensuring that solutions exist only for certain discrete values of 
A. Here however we can expect to find a solution for a range, or ranges, of A. Therefore if we have 
any eigen-operators, we have uncountably many, and operator dimensions are not quantized. 

If this is the real physics of the situation, fixed points for quantum gravity are unlike any 
fixed points discovered before in quantum field theory. Furthermore the situation is hopeless for 
asymptotic safety, since we now have an uncountable infinity of fixed points each with either no 
renormalised coupling constants or an uncountable infinity of them. 

Fortunately, sensible and familiar fixed point behaviour can be recovered if we find another 
constraint to reduce by one the dimension of the parameter space. Although the equations were 
derived by assuming a space of constant positive curvature (the four-sphere), evidently quantum 
gravity has to make sense also for negative curvature. Furthermore, one should be able to pass 
smoothly from one to the other, since a four-sphere with vanishing curvature, should be physically 
indistinguishable from flat space or a hyperbolic space with vanishing curvature. Therefore for 
both mathematical and physical reasons, we see that we need to extend the equations so as to 
incorporate consistently spaces of different topology, or at least to incorporate spaces where the 
local curvature can have either sign. 

Although the f(R) approximation so far derived has assumed the background space is a four- 
sphere, we explore the above proposal by considering its analytic extension to R < 0. Indeed in 
this case, we get one further singular point: at R = i?_ = — R + [17] . We should now thus expect 
at most a discrete set of fixed point solutions, and around these we should expect to find a discrete 
set of eigen-values and their associated eigen-perturbations. 

In sec. [7J we explain in detail our numerical techniques. The existence of some large numbers 
(in particular twice the volume of the four-sphere 2V = 768tt 2 /R 2 ~ 7580/i? 2 ), together with 
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the singular points, which control an intricate pattern of partial solutions as we uncover, makes 
the non-linear fixed point equation challenging to analyse numerically. We show how to test that 
sufficient accuracy has been achieved, and how this relates to the various numerical tolerances. As 
we will see, expanding into some regions of parameter space requires increasing accuracy; we have 
had to utilise 25 significant figures. After extensive such numerical searches however, we uncovered 
four lines of solutions in the domain i?_ < R < but none of these extend to R = — oo. Although 
at first sight, a discrete set of four solutions match onto those found for R > 0, we will see that, 
after a careful analysis of these four solutions, only three of them actually do. 

We should not be too quick to draw conclusions from the null result however. Expecting the 
equations to make sense for R < is probably too optimistic, especially given the nature of the 
approximations used to derive them. What is required is to formulate equations that correctly 
take into account throughout that R can be negative, however this task is beyond the scope of the 
present paper. 

Parameter-counting is only the beginning of a mostly-asymptotic analysis of the LPA in scalar 



field theory, which allows one to prove many of the expected features [22 23 . Perturbations 
around the fixed point either have real quantised eigen-values and power-law behaviour at large (p, 
with a power that follows from dimensional reasoning (these are the eigen-perturbations), or the 
perturbations grow faster than a power at large <p. One can prove that the latter perturbations play 
no role since one can show that any solution V(<p, t) can be written as a convergent expansion over 
the eigen-perturbations at any scale below the overall ultraviolet cutoff. One can prove that only 
the (marginally) relevant perturbations that grow as a power of the field at large ip, are associated 
with renormalised couplings. One can also prove that if the physical scale introduced by these 
renormalised couplings as they leave the fixed point, is finite, the evolution in t to the far infrared, 
tends to a finite limit (at least in the regime of large <p). In sec. [6j we consider the present case. 
Although all perturbations seem to behave in a reasonable way in the asymptotic regime, the fact 
that quantum corrections do not decouple here limits the power of these arguments. We will show 
that small finite irrelevant perturbations shrink into the fixed point, although not quite in the 
way that would naively be expected, and that ultimately the evolution under t of finite relevant 
perturbations can be expected even in the asymptotic regime to be controlled by the equations at 
finite R. 

Finally in sec. [8j we discuss some of the open questions and our findings in more detail, and 
highlight our main conclusions. 
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2 A very brief review 



Here we introduce the f(R, t) renormalisation group equations in order to set the stage for their 
analysis in subsequent sections. The purpose of the review is to expose as briefly as possible the 
key assumptions and methods that go in to deriving them. It is not our intention in the present 
paper to explore these in detail although we do make some comments here and in the conclusions. 
The starting point for the derivation of the f(R) approximation is the exact renormalisation 



group flow under the cutoff k 11,25 , written in terms of the Legendre effective action 23,26-281 



where it is commonly known as the effective average action |27j. For a single scalar field cp, it takes 
the form: 



dk 



1 



tr 



8 2 T 

6<p5(p 



dk 



71. 



(2.1) 



It is the continuum expression of Kadanoff blocking, the first step in constructing a Wilsonian 



renormalisation group 11 . The computation is performed in Euclidean signature, tr standing for 
a functional trace over the space-time coordinates. 7Z is some infrared cutoff function written as a 
momentum dependent effective mass term hip ■ 1Z ■ ip. Its purpose is to suppress momentum modes 
p < k. Evidently dependence of 7Z and thus T on k, should be understood, even though we do not 
indicate it explicitly. 

This equation needs to be significantly adapted to cope with quantum gravity, but the basic 



features, the inverted two-point function and differentiated cutoff function, remain. Therefore (2.1) 
can be interpreted schematically as the flow equation for quantum gravity, if ip is taken to stand 
for all the required fields. As we will see shortly, not only will we need the metric on some 
Euclidean-signature manifold, but also associated ghost fields and auxiliary fields. 



In the approximations we will be discussing 15-17 , the first step is to use the background 
field method, where a background field g^ u is split off by writing g^ v = g^ v + h^ u , the second 
order functional differentiation in (2.1) now being performed with respect to |2j. Gauge fixing 
and cutoff terms for h^ v are introduced in a way that leaves diffeomorphism invariance for the 
background field undisturbed. 

Although some computations are carried through in general dimensions, using various choices 
of gauge fixing and cutoff functions, we will work exclusively in four dimensions, since the final 
equations for f(R,t) are only reported for this case, and similarly we will quote only the final 



choices of gauge fixing and cutoff functions used in refs. 15-17 to derive these equations. 



S 



The gauge choice implemented is a limit such that Landau gauge with 



(2.2) 



is satisfied, where is the covariant derivative with respect to the background field. Here h = h^^, 
and indices are raised and lowered using the background field. With the gauge fixing come ghosts. 
They appear bilinearly in the action and these kinetic terms in turn need a cutoff to suppress their 
modes p < k. In standard fashion, diffeomorphism invariance of the total field g^ v becomes BRS 
invariance, which however is then broken by the cutoff. In principle it can be recovered once the 
flow k — > is complete, providing modified Ward identities are satisfied |2j. 

In order to facilitate the computation of the inverse of the two-point function, a transverse- 



traceless decomposition of the metric fluctuations is made [29 - 32 

h 



h 



1 



(2.3) 



where £ M and a are the gauge degrees of freedom to be distinguished from the physical traceless- 



transverse mode h^ v and physical scalar mode, and 



h T>1 



0. 



V^ = 0, v% 



0. 



h=h- VV . 



(2.4) 



The vector ghosts are also decomposed into their longitudinal and transverse parts. These changes 
of variables have associated Jacobians which then become bilinear kinetic terms for some (unphys- 
ical) auxiliary vector and scalar fields. These in turn need cutoffs to suppress modes p < k. 



The authors of refs. 15 - 17 then make three types of approximation. Firstly, k dependence in the 



ghost |2| and auxiliary terms on the left hand side of (2.1) is neglected. This means that the ghost 



and auxiliary field terms only contribute in (2.1) through inversion of their two-point functions 



Secondly, mixed terms depending on both and g^ u in the left hand side are neglected. This 
means that the h^ u functional derivatives in (2.1) can be regarded as evaluated at hup = 0, and 



afterwards the background field g^ u and full metric g^ u may be identified. Finally, the dependence 



on the full metric is truncated to that of (1.1 ), which means that on evaluating the right hand side 



of (2.1), only terms that take this general form are kept. 



Since it greatly helps in moulding and extracting such terms to specialise to a maximally 
symmetric space, all the authors choose to work on a sphere. In addition all authors choose to 



work with Litim's optimised cutoff profile 24 since this also leads to technical simplifications: 



r(z) = (1 - z)6(l - z) . 



(2.5) 
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In particular in the numerator in (2.1), one appearance of this can appear differentiated, thus as 
dr/dz = -6(1 - z). 



From here on there are differences in implementation. We discuss the earliest case, ref. 15 



first. It uses a cutoff of "type lb". The "b" refers to the use of (2.3), while "Type I" 16 refers to 
the fact that in (2.1), 1Z is chosen so that it modifies all appearances of the Laplacian A := — V 2 



in the two-point functions by the replacement: 

A^ A + k 2 r(A/k 2 ) 



(2.6) 



Combined with (2.5), it can be seen that the numerator restricts the eigenvalues of A to be less 
than k 2 , and in the two-point function itself this means that A is just replaced by k 2 . 



The space-time trace in (2.1 ) is taken from a heat-kernel asymptotic expansion, however certain 



unphysical modes for the component fields in (2.3) are then explicitly subtracted, namely those that 



give zero in (2.3) and thus do not contribute to These correspond to the constant mode of 



a, vectors £^ satisfying Killing's equation, and since h is taken as the physical scalar component 
mode, also the modes V^cx that satisfy the conformal Killing equation 



(2.7) 



These excluded modes correspond to zero modes for the auxiliary field two-point functions which 
are also subtracted. A similar procedure is carried out for the vector ghost fields |15| . 

Finally the rescaling step of Wilson's renormalisation group is achieved by writing everything 
in dimensionless variables using appropriate powers of the cutoff k. From here on that is what we 
will mean by / and R. If we rename the dimensionful curvature and Lagrangian by using a tilde, 
these dimensionless variables are given by 



R := R/k 2 and f(R) := f(R)/k 4 . 



(2i 
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The result is: 



384vr 2 (d t f + 4/ - 2Rf) = 
{12 + 4R 

WR 2 9(1 - %)-R 2 6(l + 



(2.9) 



5 R*e (i - f ; 



oo 



R 2 ) 



+ 



36 + 61? 



67 r> 
60 



2 



(ft/' + 2/ - 21?/") (10 - 51? - ifi? 2 + ifi? 3 ) + /' (60 - 201? 



271 

18 



R 2 



[/ + /'(! 



5ft 2 
2 



(ft/' + 2f - 2Rf") {r(-f ) + 2r(-f )} + 2/'0(l + f ) + 4/'0(l + f )] [/ + /'(l - f )] 1 

(ft/ + V' - 2Rf")f (6 + 3i? + §R 2 + Jfe I? 3 ) 
+ (ft/»-21?/''') (27-^-1^-^,^) 



+f (216 - f R 2 - ff i? 3 ) + /' (36 + 12R + §R 2 ) } hf + 3/(1 - + 9/" (1 



fl^2 
3 ' 



where 



S = 10^0 (1 - f ) . (2.10) 
Here / means f(R, t), and primes indicate differentiation with respect to R. The explicit appearance 



of r from (2.5) and the Heaviside 8 functions, arise via isolated modes, either directly or as a 



consequence of re-expressing the heat kernels in terms of unconstrained fields. The first two lines 
of the right hand side are independent of f{R) and encapsulate the contributions from the ghosts, 
auxiliaries, £^ and a. The third and fourth line arises from hT v , whilst the final ratio is the 
contribution from h. 

In ref. [l6] , the same cutoff is used but the equation is adapted to polynomial truncations only, 
which means that the Heaviside 9 functions are all set to one. Unphysical modes are isolated 
differently. The authors first consider a way of accounting for these excluded modes which amounts 
to setting £ = |32|. Then for a different way of isolating modes, they supply results for two ways 
of implementing the Landau gauge limit (the second involving higher derivatives in the 't Hooft 
averaging procedure). The resultant effect is to change S to one of three optionsj^] 

10R 2 (R 2 -20^ + 54) 10(1112-36) 



0, 



or 



(2.11) 



(R-S)(R-i) (R-3)(R-4)' 

We should make some further comments on the use of the cutoff profile ( |2.5[ ) together with 
heat-kernel methods of evaluation as in refs. |15l|16|. As noted in footnote 15 of ref. (15|, it results 



4 According to refs. 15, 16 there are not meant to be any further differences, however the equations actually differ 
by a volume factor, and also differ in the coefficient of the R 3 (d t f' + 2/' — 2Rf")/ (3/ + 3/" — Rf) term. Presumably 
these differences are typographical errors. They do not alter our conclusions. 
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in appearances of 6( n \k), which were neglected. It is not clear how these should be treated or 
interpreted. It also results in one of the major simplifications in that the heat kernel expansion is 
then truncated at a finite order. However the heat kernel expansion is an asymptotic expansion for 
small R^vaX, which means that we can expect it to have zero radius of convergence. The fact that 
the cutoff then truncates it to finite order does not therefore mean that the result is exact. In fact 
the result is smoothed implicitly by this procedure as remarked in ref. [17] . The exact result has a 
stair-case behaviour as we will shortly review. Although in principle divergences render space-time 
traces ambiguous, traces that are properly regulated by the cutoff, have an exact answer, and this 
is the one in principle that should be used. In practice using some smoothing approximation is 
presumably as legitimate as the truncations imposed in order to make progress with the exact 
renormalisation group equations in the first place. 



Benedetti and Caravelli 17 made a number of significant advances. They define Ao := A — R/3, 
Ai := A — R/A and A2 := A — -R/6, for scalar, vector and traceless-tensor component fields 
respectively. They observe that the two-point functions for the unphysical scalar and vector fields 
have factors of Aq and Ai respectively, causing them to vanish whenever the corresponding Aj has 



a vanishing eigenvalue. After replacement (2.6) and rescaling (2.8), these two-point functions thus 
have a factor of 1 — R/S (respectively 1 — i?/4). These denominator factors then appear through 



the isolated modes as poles in the first two lines of the right hand side of (2.9). 



To eliminate these singular points, they use a cutoff of type II [16] in the following sense. For 
the scalar, vector and traceless tensor component fields, the Laplacian is eliminated in favour of 
the Aj respectively. Then the cutoff 1Z is chosen differently for these component fields so that in 



their two-point functions the replacement (2.6) is made not for A but for Aj. 

The ghost action is chosen following ref. [33] . It is formally equivalent to the standard im- 
plementation but has in particular the merit of realising exact cancellations between ghost and 
gauge degrees of freedom in the Landau gauge. This greatly simplifies the contribution from the 



unphysical scalar and vector sectors. In (2.3), h is chosen as the physical scalar field in preference 
to h. This means that the a modes satisfying (2.7) now do contribute to h^. This in particular 
has the merit of removing a further singular point from the renormalisation group equations. 



The space-time traces are computed exactly by recognising that use of (2.5) on a four-sphere 
collapses this evaluation to a sum over certain powers of the eigenvalues fc 2 Aj of Aj up to the 
maximum Aj < 1 (where we have written the eigenvalues in scale- free form). However in order 
to avoid the resulting "staircase" behaviour, the results are smoothed by replacing them with the 
simplest functions matching the leading behaviour at small R and also as R — > 00. Multiplying 
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through by R 2 , the result is 



17 



384tt 2 (dtf + 4/ - 2Rf) =T 2 + T! + T^ v + 7? , (2.12) 



where the contribution from is 



from the vector modes is 



20 (d t f -2Rf" + 8f) 
72 ~ (i?-2)/'-2/ ' (2 - 13) 



71 = -36, (2.14) 



from the non-physical scalar modes: 

7^ np = -12-5i2 2 , (2.15) 

and finally from h: 

7-h 
— 

(R 4 - 5AR 2 - 54) (dtf" - 2Rf") - (R 3 + 18R 2 + 12) (dtf - 2Rf" + 2f) - 36(R 2 + 2){f + 6/") 

2(-9/" + (i?-3)/'-2/) 

(2.16) 

Again, / means f(R,t), and primes indicate differentiation with respect to R. 

Although the singular behaviour at R = 3 and R = 4 has thus been avoided, singular behaviour 
still remains in the h sector although for different reasons. These singular points appear as algebraic 
zeroes in front /'". 



Firstly, introducing the cutoff 1Z in order to modify the Laplacian as in (2.6), using either A 



or the Aj, means that 1Z is forced to depend on the other terms in the two-point function. In the 



h sector this includes an /" term. In (2.1), when the cutoff gets differentiated with respect to k, 



and after the change of variables (2.8), this turns into the dtf" — 2Rf" term present in (2.16) and 
the penultimate line of (2.9). As Benedetti and Caravelli have noted [17], the modification (2.6) 
is thus responsible for turning a second-order differential equation into a third-order one (through 
introducing /'"), and at the same time removing the extra parameter by introducing a singular 
point at R = (through the — 2R coefficient). 

Secondly, the lowest eigenvalue in the h sector corresponds to the constant mode. For the 
modified Laplacian Ao, this eigenvalue is now negative: Ao = —R/3. This mode is then never 
eliminated by the cutoff requirement Ao < 1. The f"(R) piece turns out to be multiplied by 

Ea,,<i( 1 ~ A o)- Note that the 

negative eigenvalue contributes increasingly negatively to this sum 
13 



once 22 > 3, whilst the contribution of the positive eigenvalues falls to zero as 22 increases. Thus 
there is always a positive value of R for which this sum changes sign. After smoothing, this turns 
in to the factor 22 4 — 5422 2 — 54 in (2.16), which changes sign by vanishing at the points R = R±, 
where 

R+ = -R- = \J27 + 3V87 . (2.17) 



3 Basic requirements on the solutions to the fixed point equations 

To obtain the corresponding fixed point equations we simply set dtf = 0, and / now stands for 
the ^independent fixed point solution f(R). As already advertised, these are then third order 
non- linear ordinary differential equations, the first being 



768tt 2 (2/ - Rf) = 

5R 2 6(l- f) - (12 + 422 - § 22 2 ) 



(3.1) 



+ 
+ 
+ 
+ 



10 R 2 9(1 - f ) - R 2 9(1 + f ) - (36 + 6 22 



67 r>2\ 
60 K ) 



1-f 



-1 



(2f - 2Rf) (10 -5R- ?§R 2 + gf 22 3 ) + /' (60 - 2022 - ?§R 2 ) ] [/ + f'(l 
(2f - 222/") {r(-f ) + 2r(-f )} + 2f'6(l + f ) + 4/0(1 + f )] [/ + f'(l - 
(2f - 222/")/' (6 + 322 + f 2? 2 + ^ 22 3 ) 

22 4 ) 



3 > 



-2Rf" (27 



91^2 29 e)3 



R" 



181 D 4^ 
3360 



+/" (216 - f 22 2 - §22 3 ) + /' (36 + 1222 + |22 2 ) 1 hf + 3/'(l - §22) + 9/"(l 



R\2 
3 > 



where S is given by (2.10). The second is given by the above with all 6 functions set to one, and 



S given by (2.11). And the third is given by 

768vr 2 (2/ - 22/') = T 2 + Ti + T* p + 7? , 



(3.2) 



where 



To 



40 (22/" -4/Q 
(22 -2)/' -2/' 



and 



~- h _R(R 4 - 5422 2 - 54) /'" - (22 3 + 1822 2 + 12) (22/" - f) + 18 (22 2 + 2) (/' + 6/") 



T' 
'n 



(3.3) 



(3.4) 



9f" - (22 - 3)/' + 2/ 

Since the flow equations were derived on a four-sphere it is natural to assume that the domain 
over which solutions should exist is < 22 < oo, where we include 22 = via the limit 22 — > 0. 
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As already mentioned in the introduction, we should insist that solutions to (3.1) or (3.2), and 
all derivatives of the solutions, exist throughout this domain. Since these fixed point solutions 
already give us massless continuum limits, failure to satisfy these requirements would translate into 
unacceptable singularities in the physics implied by these fixed points. Also these singularities would 
be inherited by the more general field theories which would need to be built using the renormalised 
trajectories emanating from these fixed points]^] Note that the limit R — > is included in the 
domain, because a fundamental requirement of the Wilsonian renormalisation group is that the 
effective Lagrangian f(R,t) must be quasi-local (for non- vanishing k). This means that for non- 
vanishing k, the effective Lagrangian must have a derivative expansion to all orders. Since R 
contains two derivatives, this requirement implies that f(R,t) must have a Taylor expansion in R 
about R = 0. 

Again, as touched on in the introduction (and see also the conclusions), these requirements can 
fail to be satisfied, not because of underlying difficulties in the physics, but because a poor choice 
of cutoff results in an ill-defined Wilsonian renormalisation group. Generally in this case it will not 
be clear how to match across the singularities caused by the cutoff; we then have no option but to 
abandon this choice of cutoff and search for a better one. 



The first fixed point equation, (3.1), is a case in point. The appearance of explicit Heaviside 



9 functions means that any fixed point solutions cannot possibly have well defined derivatives to 



all orders. Clearly this is inherited from the fact that the cutoff profile (2.5) itself incorporates the 
8 function, resulting in it having singular second and higher derivatives. However, we do not yet 
have to abandon this equation. We can continue if we assume that the correct matching across 
a 9(R — R c ) is provided by absorbing this change in a jump in f"'{R) at R = R c . The equation 
can then be solved separately in the two regions R < R c and R > R c , thus also leading to discrete 
changes in all the higher derivatives as we cross the boundary from one region to the other. 

4 Counting parameters in the fixed point solutions 

We are now ready to address a central theme of the paper: the counting of parameters in the 



solutions to (3.1) or (3.2). We start by recalling some textbook mathematics adapted for the 



present purpose. We will often need to consider these equations cast in normal form: 

f'"(R)=rhs, (4.1) 



3 See for example further discussion and references in the review ref. 23 about all these issues 
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where rhs contains no /"' terms. A Taylor expansion about some generic point R„ takes the form: 



f(R) = f(R p ) + (R- R p )f'(R p ) + \(R- R P ) 2 f"(Rp) + l( R - Rpff±(R P ) + ■■■, 



(4.2) 



where the differing values for f"'{R p ) and higher derivatives are required if R p coincides with the 



jump in a 9 function. Since (4.1) determines the fourth coefficient in terms of the first three, we 



see that typically (4.2) provides a series solution depending on three continuous real parameters, 
here 

f(Rp), f'(Rp) and f"(R p ), (4.3) 

with some finite radius of convergence p whose value also depends on these parameters. Therefore we 
recover the standard result that around a generic point R p there is some domain T> = (R p —p, R p +p) 
in which we have a three-parameter set of well-defined solutions. From here we can try to extend 
the solution to a larger domain, e.g. by matching to a Taylor expansion about another point within 



T>. A typical problem, seen also in the LPA and the derivative expansion 23 , is that eventually. 



at some point R = R c , dependent on the parameters, the denominator of rhs develops a zero, so 



that as R — )• R c , (4.1 ) implies 



f"(R)~2c ± /(R-R c 



(4.4) 



where is also some function of the three parameters (again two values are required if R c coincides 
with a jump in a 9 function) and "~" means we are neglecting less singular parts j^] Thus we see 
that the solution typically ends in a moveable singularity, of form 



f(R) ~ c ± {R - R c ) 2 In \R - R c \ . 



(4.5) 



The fixed point equations (3.1) and (3.2) present another challenge in that they also have fixed 



singular points. These correspond in rhs to explicit algebraic poles in R. 



4.1 Counting parameters of the fixed point solutions of refs. 1 15 , 16 



Looking at the structure of (3.1 ), we see that from the first two lines there are single poles at R = 3 



and R = 4, as already discussed in the previous section. It is straightforward to check that they 
are not cancelled by factors in the numerator on either side of the relevant 9 functions, i.e. the 
first line does indeed diverge for R — > 3 from above, and also from below, and similarly the second 
line diverges for R — > 4 =t . /'" appears once in the penultimate line. Rearranging the equation into 



3 Exceptionally, under appropriate conditions higher powers of the linear factor R — R c can appear in (4.4 1 
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normal form it is clear by inspection that the explicit single poles at R = 3, 4 still do not disappear, 
and indeed we also gain poles from inverting the polynomial 

R („-nje-»#-™#), (4.6) 

V 20 30 3360 / ' y ' 

since this appears as a factor multiplying /'". The term in brackets has four single zeroes and one 
of them is real and positive: R = 2.0065. We therefore have four fixed single poles in rhs in the 
applicable domain, namely at R = R c = 0, 2.0065, 3, 4. 



At these points, / will end at a singularity of form (4.5) unless the /-dependent parts in rhs 
conspire to cancel the pole. Substituting the Taylor expansion ( |4.2[ ), with R p = R c , we see that this 
requirement forces some generally non-linear combination of f(R c ), f'(R c ) and f"(R c ) to vanish. 
Actually, at two of the singular points the constraints simplify to become linear: 

/'(3) = |/(3) and /" (4) = 5/' (4) - 2/(4) , (4.7) 

but the constraint at the origin is: 

2 192^ 2 /(0) 3 + 6/(0) 2 + 480vr 2 /(0) 2 /'(0) + 2/(0)/'(0) + 288^ 2 /'(0) 2 /(0) - 9/'(0) 2 



/"(0) 



9 3/(0) + 192tt 2 /(0) 2 + 192tt 2 /(0)/'(0) - 7/'(0) 

(4- 



and the constraint at R = 2.0065 is also non-linear and even longer. As we argue in the next 
subsection, these four constraints provide four independent boundary conditions on the three pa- 
rameters in our solution. Therefore the equation is over-constrained and it is not credible that any 
non-singular solutions exist over the range < R < 4 let alone over the full domain < R < oo. 



Choosing £ to be one of the expressions in (2.11), setting all 9 functions to one, and repeating 



the analysis, we still have the same number and positions for the fixed singularities and moreover 
these yield exactly the same constraints. Therefore the second equation also is over constrained 
and thus has no smooth solutions in the range < R < 4. 

4.2 Why the counting arguments should be trusted 



If the equation (3.1) could nevertheless furnish a smooth solution compatible with these four con- 
straints then it must be that somewhere in parameter space these constraints are no longer inde- 
pendent. Consider two fixed singularities bounding a domain T> which is free of fixed singularities 
(e.<7- 3 < R < 4). If the corresponding constraints are not always independent, instead of obtaining 
the expected one-parameter sets of solutions within T>, we will also find regions of parameter space 
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where there are still two-parameter sets of solutions. We now explain why this is too much to ex- 
pect from these equations. It is straightforward to extend the arguments to the full set of singular 
points and this is one reason why we are confident in these counting arguments, and the conclusion 



that (3.1) has no global solutions. 

From the previous subsection we know that in the neighbourhood of each of the fixed singular 
points bounding D, we have an analytic solution (viz. a convergent Taylor expansion) which, in 
order to avoid a singularity in f"'(R), has only two parameters. Let us choose a generic point R p 
inside T>. If we are lucky, we can extend our two-parameter solution at the upper singular point 



down to R p . In the space of parameters at R p , namely (4.3), we see that the set of such solutions 



form a two-dimensional surface. If we are unlucky, the solution will fail at a moveable singularity of 



form (4.5 ) before we reach this point. Therefore this two-dimensional surface has boundaries (which 
may thus also result in disconnected pieces). Now we develop the solution up to R p from the lower 
singular point. This will result in another two-dimensional surface with boundaries. Given that 
there is nothing special about the point R p , unless there is some hidden simplicity or symmetry 
in the equations^ these two surfaces, if they meet at all, do so on surfaces of one less dimension, 
i.e. on curves (with possibly disconnected pieces). Note that we can place the point R p anywhere 
inside the domain, but the topology of the space of matching solutions must be independent of 
this choice. The only way a solution space of dimension two could then result, is if the equations 
already have this regularity built in across the whole domain. Although we will not directly need 
it, we note that these remarks apply equally well if one of the singular points is sent to R = oo. 

Let us explain further by considering a counter-example. Consider the third order differential 
equation: 

r = R2f " ~ 2Rf + 2/ (4 9 ) 

1 R(R-l)D(f,f>,f»,R)> { -> 
where D is any combination of its arguments assumed to be well behaved in the region < R < 1. 
This region is bounded by two fixed singularities which impose on a smooth solution the conditions: 

/(0) = and /"(l)-2/'(l) + 2/(l) = 0. (4.10) 

Naively such a solution has thus only one parameter in this domain, however the solution actually 
has two parameters (a, b) and is given by the vanishing of the numerator in rhs throughout the 
domain: 

f{R) = aR 2 + bR. (4.11) 



7 hidden because such regularity is clearly not evident in (3.1 1 



18 



Of course the reason this is allowed as a solution to (4.9) is that this is also separately a solution 
of the left hand side: f"'(R) = 0. This is an example of the type of hidden simplicity that would 
be required to evade the counting argument. 



We can check that this sort of factorisation does not happen in (3.1), for example functions 
given by the vanishing of the final denominator over some range of R: 

2f + (3-2R)f + (3- R) 2 f" = 0, (4.12) 

do not also satisfy the final numerator. If this were the case, these functions would be solutions of 



the normal form (4.1). Of course, there is nothing in the derivation [15, 16 which would lead one 
to expect such a remarkable coincidence, and had it been the case it would point to a much deeper 
understanding of the renormalisation group in quantum gravity than currently gained from these 
fixed point equations. 

In the next subsection, we will see what to expect from the latest fixed point equation ( |3.2[ ) |17| . 
We will then follow this up with an extensive numerical investigation in sec. [7] from which we will 
confirm these expectations. This is another reason for our confidence in these counting arguments. 
We will also show how easy it is to be misled by the numerics if one is willing to abandon these 
counting arguments too quickly. 

4.3 Counting parameters in the fixed point solutions of ref. |17| 



The fixed point equation is (3.2), where the terms on the right hand side are given by (3.3), (2.14) 
(2.15) and (3.4). /'" appears once: in (3.4). After casting in normal form (4.1), and using (2.17) 



we see that we have just two fixed singular points corresponding to the algebraic single poles at 
R = and R = R+. Since we take the domain of the solution to be < R < oo, the algebraic 
single pole at R = R- need not concern us. Requiring that a Taylor expansion at R = exists, we 
find that /"(0) is determined in terms of /(0) and /'(0) as: 

2 6a 2 - 9ai 2 + 4807r 2 a 2 ai + 2a ai + 1927r 2 a 3 + 288vr 2 a ai 2 
a2 ~~9 3a - la\ + 1927r 2 a 2 + 192vr 2 o ai ' ^' ' 

where it will be convenient to define a n := f^ n \0). Similarly requiring a Taylor expansion at 

R = turns out to imply a quadratic constraint on f"(R + ) in terms of f'(R+) and f(R+). As 

before, we expect these singular points to provide independent boundary conditions, therefore we 

are left with at most a one-parameter set of solutions - with possibly disconnected pieces. 

As we discussed in the Introduction, another constraint is usually supplied by the asymptotic 

behaviour at large field. As we will see from sec. [5] however, the asymptotic behaviour does not 
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result in a reduction of parameters. Instead we find that as R — > oo, 

f(R) = AR 2 + rS^A + Bcos\nR 2 + CsmlnR 2 

o 1 728 1 92 

- — (103B - U9C) cos In R 2 - J—-A(11C + 7B)tt 2 cos In R 2 - — (B 2 + Y1BC - C 2 )tt 2 cos In R 4 

o i 798 384 

- — (1495 + 103C) sin In R 2 + -^-A{-7C + HS)vr 2 sin In R 2 + — (35 2 - BC - 3C 2 )vr 2 sin In i? 4 

+ —A - 96(9^ 2 + 2B 2 + 2C 2 )vr 2 + 0(1/ R) , (4.14) 

where A,B,C are three real parameters subject only to the constraint that they lie within a cone 
given by the inequality ( 5.17[ ) (otherwise f(R) develops periodic singularities, as shown in sec. [5]). 
Therefore we should still expect one-parameter set(s) of solutions, i.e. lines of fixed points. As we 
will see in sec. [7[ that is precisely what we find. 

If true, this would be a surprising outcome for the asymptotic safety scenario. We would need 
to use other arguments to decide which line of fixed points corresponds to the one realised by 
quantum gravity. On the chosen line of fixed points, we would have an exactly marginal coupling 
(i.e. dimensionless non-running parameter) labelling the position of the experimentally realised 
fixed point on this line, and an exactly marginal operator corresponding to perturbations that 
move the theory along this line. We would need to use experimental measurements to determine 
the value of this coupling. The real problem with this conclusion however is that, as we will see 
in sec. [6| it implies that the renormalisation group eigenspectrum is not quantised, i.e. each fixed 
point supports a continuous infinity of operators with a continuum of scaling dimensions. 

As already argued in the introduction, the resolution must lie in the fact that we have not 
yet considered the full domain over which quantum gravity must be defined. Although the flow 
equations in sec. [2] were derived by working on a four-sphere, quantum gravity must also make 
sense for spaces where R < 0. In fact our insistence on quasi- locality (cf. sec. [3]) means that our 
solution f(R) already exists in some small neighbourhood of R = 0, and thus also for R < to 
some extent. As explained in the introduction, we will therefore explore this option in this paper 



by analytically continuing the flow equation (2.12) to R < 0, i.e. we will regard this equation and 



the corresponding fixed point equation (3.2), as valid for all real values — oo < R < +oo, despite 
the explicit use of a four-sphere in their derivation. 

It is easy to see what happens to the parameter-counting arguments now. The singular point 



at R = R-, cf. (2.17), imposes a quadratic constraint on f"(R_), determining it to be one of the 
roots in terms of f'(R-) and This constraint is independent of the one at R+ because the 

fixed point equation (|3.2l) is not symmetric under R i— >■ —R. The three singular points R = 0, R± 
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together thus ensure that we have at most a discrete set of solutions. These solutions are however 



subject to the asymptotic constraint (5.17) which now operates independently on both sides, and 



each branch can be expected to eliminate some of these solutions. We have checked that the 



asymptotic expansion (4.14) is valid also for R < 0. However, because (3.2) is not symmetric under 
R h-> —R, there is no simple relation expected between the values of the asymptotic parameters 
A±, B±,C± for positive/negative R, cf. sec. [5} 

Although the asymptotic behaviour is unusual, we have arrived at the expected conclusion that 
there exists at most a discrete set of global solutions. We can hope that, when properly derived and 
understood, the equations furnish just one solution: the fixed point about which asymptotically 
safe quantum gravity can be constructed. As we will see in sec. |6j parameter-counting now shows 
that the eigen-spectrum behaves properly too: the renormalisation group eigenvalues are quantised. 
Therefore we can also hope that the relevant ones are only finite in number. 

Unfortunately, as described in sec. [7j despite an extensive numerical investigation we found 
only partial solutions and no solution defined for all R. 

5 Asymptotic expansion of the fixed point solutions 

We have seen in sec. [3] that the solution must be Taylor expandable about R = as a consequence 
of the requirement of quasi- locality. It is for a very different reason that constraints arise as R — > oo. 
We have already noted that partial solutions typically end at a moveable singularity, this being of 



forrrrl (4.5). This is the typical situation one finds in the LPA and derivative expansions 23 . The 



mere fact that the solution is well defined for all R > Rf/s, where Rf is some fixed finite value 
and e is small, would therefore be expected to impose some strong constraints on the solution. 
Furthermore, one can hope that now we are supplied with the small parameter e, we can develop 
the solution analytically. This is the asymptotic expansion for f(R). 

Note that there is no reason to expect that f(R) has a Taylor or Laurent expansion in 1/R as 
R — )■ oo. Indeed in the LPA for a single scalar field (p in two dimensions the result is a sinusoidal 
V(tp) corresponding to critical Sine-Gordon models |21], whilst at higher orders in the derivative 



expansion, V(<p) for large <p grows as some irrational power of <p 20,21,23 . In both these cases, 
V(<p) has an essential singularity in 1/ip, at ip = oo. 
8 the ± labels are not applicable now 
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Benedetti and Caravelli found the asymptotic expansion 17 : 



f(R) = AR 2 + *AR +^A- 864vr 2 A 2 + 0(1/ R) , (5.1) 

(and computed to high order further coefficients as a Taylor series in l/i?)j^] From this it would 
appear that requiring a fixed point solution to exist for large R results in two conditions on the 
third order differential equation (and thus only one parameter A). 

So far this is very similar to the behaviour found in the LPA in greater than two dimensions: 
the leading power being determined by dimensions as would be expected. However there is still 
an important test that should be carried out [19| . By linearising around the solution, one can 
understand what has happened to the 'missing' parameters. In the LPA for scalar field theory, 
one finds that the asymptotic expansion indeed describes an isolated one-parameter set because 
perturbations corresponding to the missing parameter, grow much faster than the asymptotic ex- 
pansion (in fact as the exponential of a power of cp). No analytic or numerical solution is found 
that corresponds to such behaviour at the non-linear level. This can be understood heuristically 
because the linearised solution for such behaviour involves balancing the same terms that result in 
the moveable singularities, and thus can be expected to suffer this fate at the full non-linear level. 

Carrying out this test on the present case, we will arrive at very different conclusions. Writing 



f(R) i — y f{R) + fif(R) m (3-2 ) and expanding to first order in 5f gives us just the right hand side of 



( |6.1[ ) (with the left hand side set to zero). Substituting the expression (5.1) for f(R) (actually we 
need only the first two terms) and keeping only the leading terms in large R in front of the 5f( n \ 
we obtain the linear ordinary differential equation: 

R 3 Sf" - R 2 Sf" + 6R Sf - 10 Sf = . (5.2) 

This equation is invariant under scale changes R t— >■ sR and thus has power law solutions. The 
complete solution is given in terms of a linear combination of powers R p , where p = 2, 1 + 2i, 1 — 2i. 
Expressing the complex powers in terms of their real and imaginary parts: 

5f(R) = 5AR 2 + 5B R cos In R 2 + 8C Rsin In R 2 , (5.3) 

where we have introduced two infinitesimal parameters SB and SC, and recognised that the first 



term results from perturbing A in (5.1). We have therefore discovered that (5.1) is not an isolated 



9 Actually their numerical coefficients are wrong as a result of the error in (2.161, but their series has the same 
qualitative behaviour as the one displayed here. 
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one-parameter set: perturbations in the missing parameters correspond to sub- leading corrections 
in the large R limit. 

We also see that the full asymptotic expansion including these terms cannot be made to corre- 



spond to a Taylor expansion in 1/R. Given the complexity of (3.2) and the anticipated complexity 



of the series solution following (5.3), it will prove fruitful to separate it into the various orders of 



the asymptotic expansion by introducing a small parameter e. For this purpose define: 

f e (R) = e 2 f(R/e). (5.4) 



where R is kept finite, and we are taking the hint from (5.1) and (5.3), that lim e ^.o fs(R) = fo(R) 
exists. 

(In principle there could also be asymptotic solutions at large R other than the assumed / oc R 2 
however no other singularity-free behaviour has been found numerically. It is straightforward to 
rule out analytically various ansatze: for example suppose that at large R, f oc RP for some p ^ 2. 
If p < 2 then the R 2 asymptotic will dominate almost always (i.e. unless A = 0). If p > 2 then 



since the right hand side of (3.2) is homogeneous in /, it is easy to see that to leading order the 
left hand side must be satisfied on its own, unless the denominator of 72 or 7q vanishes to leading 
order. However satisfying any of these conditions to leading order requires p = 2, thus ruling out 
this ansatz.) 

Now we write: 

f £ (R) = fo(R) + sfr(R) + e 2 f 2 (R) + ■■■ . (5.5) 

However, because a Taylor expansion in e does not exist, the higher f n (R) actually still depend on 
e, although we will see that ultimately this implicit dependence can be absorbed in the parameters. 
The division into orders as above will still make sense because we will insist that the implicit 
dependence of f n {R) on e is not sufficient to cause it to vanish as fast as e, or to cause it to diverge 
as fast as 1/e. (This is in effect the Poincare definition of an asymptotic series.) It will be helpful 
to write 

f n (R) = R 2 - n g n (R) ■ (5.6) 



By (5.4) and (5.5), the g n (R) are then invariant under the rescaling R i— > sR and e i— > se. Note 
that the higher g n will achieve this through their implicit dependence on e. It follows that a solution 
for g n is at the correct order in the asymptotic expansion if its behaviour at large R is bracketed 
such that 

l/R<g n {R)<R- (5.7) 
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At lowest order we have go(R) = A, i.e. 



MR) = AR 2 



(5.8) 



Rewriting the fixed point equation as a differential equation for f £ {R) using (5.4), amounts to the 
substitutions R h-> R/e and /( n ) e n ~ 2 / £ (n) in jOJ, ^13| >, plo] and (jjU). We substitute 



(5.5) and expand the differential equation as a series in the explicit dependences on e. 



The lowest powers come from 0{e 2 ) terms in ( |2.15 1 and (3.4). Requiring that the 0(e 2 ) part 
vanishes gives us a differential equation for g\{R): 



R 3 g'l' + 2R 2 g'{ + \Rg' x - 4 5l = -6A . 



(5.9) 



The equation has the expected invariance under scale changes R i— > sR (explicit dependence on e 
having dropped out). After shifting g\ \-± 3 Aj '2 + g\, this equation is also solved by power ansatz, 



this time with p = l,±2i. Recalling the extra power of R from (5.6), we see that the first power 



just perturbs the leading order solution and is rejected by the conditions (|5.7l). Expressing i? ±2i in 



terms of real and imaginary parts, we recognise that we have recovered in full form the linearised 



sub-leading corrections in (5.3), thus the next order term in the asymptotic series is: 

fi(R) = Rl^A + BcoslnR 2 + CsmlnR 2 



(5.10) 



We prefer to write the sinusoidal parts in this way but note that there are a number of alternatives, 
for example we can exchange B and C for parameters E and A: 



B cos In R 2 + C sin In R 2 = E cos ln(R 2 /A) 



(5.11) 



This form makes it clear that a dimensionful parameter A = /c 4 A has crept in, which thus is 
multiplicative in the implicit e dependence and scaling (A i— > s 2 A). 



Had we not assumed the form (5.8) for the lowest order, we would have found that the dominant 
terms come from the same 0(e~ 2 ) terms as above, namely from the —5R 2 term in (2.15), and the 



terms at this order from (3.4): 



R 



. R 3 f'"-R(Rf"-f') 
2/ - Rf> 



(5.12) 



but also including the left hand side of (3.2). This limit is actually singular: substituting f = fo 



AR 2 results in the left hand side vanishing, whilst (5.12) gives 0/0. Adding in the ef\ dependence 



also, results in equations that are no longer singular, but now the left hand side of (3.2) is neglected 
because it is of higher order in e, whilst (3.4) continues to contribute at 0{e~ 2 ). Rearranged this 
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is the equation (5.9) when written in terms of g\. Although the leading part (5.8) is as expected 
by dimensions, we see that in fact the reason for its appearance is more subtle than in the LPA 
for scalar field theory. In the LPA, the reason is that the quantum corrections can be neglected 
for large field. In contrast here, the dependence at large R is actually dominated by balancing the 
physical scalar quantum corrections against the non-physical scalar quantum corrections. 

The O(l) contribution comes from the next power in e: the 0(e _1 ) terms. These are found in 
the left hand side of (3.2), and in (3.4). Changing variable to \R\ = expu, we obtain: 

{dl - 4d 2 + 9d u - 10} g 2 (u) = -315A/2 + 960(9,4 2 + 25 2 + 20 2 )vr 2 

- 3(2C + 315) cos 2u + 6912,4(5 - 2C)vr 2 cos 2u - 384(5 + 3C) (35 - C)vr 2 cos \u 

+ 3(25 - 31C) sin2u + 6912,4(25 + C)vr 2 sin2ii + 768(25 + C)(5 - 2C)vr 2 sin4u . (5.13) 

We see that the left hand side is translation invariant under u i— >■ u + In s, as required. The general 



solution is just given by (5.3), and does not satisfy (5.7). 

At higher orders n > 3 also, the linear homogeneous part is always of this form when expressed 
in terms of f n (u), and thus the general solution is always given by the perturbation of parameters 



in (5.3). It differs for the g n {u) only through the change of variables (5.6); these are all excluded 



by (5.7). It will therefore be the special solution we are interested in from now on. 



The special solution of the above equation is lengthy but readily found. Expressed in terms of 



R: 



f 2 (R) = 52(5) 
3 



63 



A - 96(9A 2 + 25 2 + 2C 2 )tt 2 



(1035-149C) cos In R l 



1728 



68 
3 

— (1495 + 103C) sin In 5 2 + 

68 v ; 17 



^(11C + 75)tt 2 cos ln5 2 



192 
"3T 



(5 2 + 125C - C 2 )tt 2 cos In R A 



17 

1728 A (-7C+ 1 15)vr 2 sin In R 2 + ^ (35 2 -BC- 3C 2 )vr 2 sin In R 4 

(5.14) 



This 0(1) part of the asymptotic expansion is the last piece that can be derived in closed form. 



The asymptotic expansion is retrieved by substituting back into (5.5) and setting the formal 
parameter e = 1. We see that, together with (5.8) and (5.10), we have now derived the already 
advertised form ( |4.14 ). 



The 0(1/5) term comes from 0(e°) in the fixed point equation. All parts of (3.2) now con- 
tribute. This gives us a differential equation of the form: 



{dl - 78 2 u + 20d u - 24} g 3 (u) = M(u)/d(u) 



(5.15) 



the modulus is included so that the derivation will be applicable also when R < 
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where the numerator M is too long to display. It is quartic in the parameters, and linear in sines 
and cosines with arguments 2nu (n = 0, 1, 2, 3, 4). In terms of R, the denominator is given by: 



d(R) = 11A/2 + (5 - 2C) cos In R 2 + (2B + C) sin In R? 



(5.16) 



The general solution of (5.15) just reproduces (5.3) for fs(R) and is disallowed by (5.7), as already 



mentioned. The special solution can be expressed in integral form. The important feature is that 
it manifests moveable singularities (4.5) whenever d(R) has a zeroj^] It is straightforward to see 
that d(R) has real zeroes if and only if the following is violated: 

121 
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■A 2 > B 2 + C 2 , 



(5.17) 



i.e. the inequality is violated if the parameters lie outside the cone with axis A, vertex at the origin, 
and opening angle 2tan _1 f^TgJ ~ 135.8°. We see that if the inequality is violated the asymptotic 
series predicts an infinite series of such singularities R c with multiplicative periodicity R c \— > R c e n , 
although only those with large enough R c can be trusted. Indeed we find that, integrating out 
numerically from the fixed singular point R + , if we can reach large enough R, we can extract 



reliable values for A, B, C by using the first three orders of the asymptotic expansion (4.14). 



At 0(1/ R n ) with n > 2, the equations take a similar form to (5.15) and can also be solved 
in terms of integrals of ratios of sine and cosine series. The denominator has a factor which is an 



increasing power of d(R), however a new factor also arises, coming from the denominator in (3.3): 

15,4/2 + (B- 2C) cos In R 2 + (2B + C) sin In R 2 . (5.18) 



Note that this is of the same form as (5.16) but is a cone with larger opening angle. 



The asymptotic series solution as R 



-co is the same i.e. as given by (4.14), developing 



moveable singularities (4.5) at the zeroes of (5.16), if (5.17) is violated. The derivation is exactly 



as above, since we have been careful to ensure that all the intermediate steps are valid also for 



R < 0. As mentioned at the end of sec. 4.3, we expect for any given fixed point solution that 
the values of A,B and C are different in the two sectors R —> ±co. As a function of complex R, 
this can be expected to arise because the asymptotic series is multivalued and thus subject to the 
Stokes phenomenon. 



1 We have verified that Af does not contain d(u) as a factor. (In fact d(u) first appears in the derivation of the #2 
equation where however it is a factor that appears also in the numerator.) 



26 



6 Analysis of eigen-operators 



In this section we develop a similar analysis to sees. [4] and [5] for the perturbations around a 



fixed point solution f(R) of (3.2). We first take the perturbations to be strictly infinitesimal and 
then go on to investigate asymptotic properties of finite perturbations. To begin, we substitute 



f(R, t) = f(R) + Sf(R, t) into (2.12) and expand to first order in the infinitesimal 6f(R, t). Writing 
Sf = dt Sf, the result is: 



384^ 2 Sf + 



20 Sf 



(R - 2)f - 2/ 



+ 



(it! 4 - 54i? 2 - 54) Sf" - {R 3 + 18i? 2 + 12) Sf 



768tt 2 (RSf'-2Sf) + 



2 (9f» -[R- 3]f + 2/) 
40 (R Sf" - 4 Sf) - f 2 ([R - 2] Sf - 2 Sf) 



(R - 2)f - 2/ 

(# 4 - 54i? 2 - 54)RSf" - (R 3 + 18i? 2 + 12) (RSf - Sf) + 18(R 2 + 2) (Sf + 6 5/") 



9/"-[fl-3]/' + 2/ 



- h 9Sf-[R-3]Sf + 2Sf 



(6.1) 



9f"-[R-3]f + 2f 

where we have utilised the expressions ( 3.3|3.4 ). Note that we have split the equation so that the 



t-derivative terms appear exclusively on the left hand side. The left hand side pieces of (2.12) give 



the first pieces appearing on the left and right of this equation, then in order we have the parts 



from perturbing (2.13) and (2.16). 



We search for perturbations with a characteristic scaling dimension under the exact renormali- 
sation group i.e. for separable solutions of the form: 



6f(R, t) = v(R) a exp (-2At) , 



(6.2) 



where here a is some dimensionless normalisation constant (and the factor of —2 has been included 
for convenience). The expectation is that the v(R) are the eigen-operators for this fixed point, 
i.e. the operators with well defined scaling dimension 12 It is also expected that the solutions 



are quantised so that we obtain a countable infinity of such operators, and that these span the 
space of all possible interactions at this level of truncation. The associated couplings follow from 



the t-dependent coefficients after untieing the transformation (2.8). Recalling that t = ln(fe//i), 
where \x is some finite physical scale, we see that they have definite mass-dimension at the fixed 
point and are given by g = a// 2A . The couplings with A < 0, A = and A > are respectively 
irrelevant, marginal and relevant. It is expected that the last, together with any marginally relevant 



12 Actually since these are perturbations of < |1.1[ ) they only yield the integrated operators. Operators that can be 
written as a total covariant derivative are excluded. 
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couplings, evolve into the running renormalised couplings which parametrise the general continuum 
limit associated to this fixed point, i.e. the most general renormalised trajectory out of this fixed 



point, and that these couplings are finite if the corresponding afi 2 ^ are finite 22 23 . Although 
these are natural expectations for a Wilsonian renormalisation group [IT] , they are not necessarily 
true in the present context. In this section we will see how far we can go in verifying or refuting 
these expectations. 



Substituting (6.2) into (6.1) we obtain: 



x . „„„ 2 40?/ (R 4 - 54i? 2 - 54V' - (R 3 + 18i? 2 + 12V 

A ^ 768tt 2 ?; + — — — - — — + v 



(R - 2)f - 2/ 9/" -[R- 3]f + 2/ 

(R 4 - 5AR 2 - 54) Rv'" - (R 3 + 18i? 2 + 12) (Rv" - v') + 18(R 2 + 2){v' + 6v") 



9f"-[R-3]f' + 2f 



- h 9v"-[R-3]v' + 2v 
9/"- [R-3]f' + 2f { ' 



This is a linear third order ordinary differential equation for v. Thus in the neighbourhood of a 
generic point R p , we can construct a Taylor expansion solution with three parameters. However, 
since the equation is linear, one of these parameters is trivial, being just the overall scale. Choosing 
some convenient normalisation, we are left with a two-dimensional parameter space. Since f(R) is 



now a fixed global solution to the fixed point equation (3.2) and v(R) does not appear in the de- 
nominators or in the coefficient of v"' , there is no opportunity for moveable singularities. Therefore 
the constraints on finding a global smooth solution for v can only arise from the fixed singularities 
and its behaviour as R — > oo. 

6.1 Asymptotic analysis of eigen-operators 



We can analyse what constraints arise as R — > oo, by developing the asymptotic solutions of (6.3). 
To find the asymptotic expansion we follow a similar analysis to that of sec. [5| We can use the 
freedom to choose a normalisation to eliminate any power of e in front of v. Therefore we define: 

v £ (R) = v(R/e) , (6.4) 

where R is kept finite, and expand as 

v £ (R) = v (R) + ev x {R) + e 2 v 2 {R) + ■■■ , (6.5) 
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where the v n also have implicit dependence on e. The combinations v n R n will be invariant under 
R i — y sR and e i— > s e, after multiplication by some power of s to be determined. [This is the missing 



power in (6.4).] We rewrite the eigen-operator equation (6.3) as a differential equation for v e (R) 
with coefficients dependent on f £ (R). This amounts to the substitutions R i-> R/e, i— > e™ -2 /^ 
and i y e n vi n \ We substitute (5.5) together with ( 5.8|5.10 5.14), and (6.5), and then expand 
the differential equation as a series in the explicit dependences on e. As before the equations are 
most easily solved by changing variables as \R\ = expu. Here we report the results exclusively in 
terms of R however. 

The lowest power comes from 0(e~ 1 ) pieces in the third terms on the left and right of (6.3), 



i.e. originally from (2.16), yielding: 



R 3 v'o - R 2 v'q + 6Rv' Q - Wv + A {R 2 v'q - Rv' } = . 



(6.6) 



As expected, this equation is invariant under scaling R \— > sR. It is therefore solved by a power 
ansatz. We find that for generic A, the three linearly independent solutions are v$ = R 2 , and 
vq = RP , where p can take two values: 

A 



P 



1 _-±IVA 2 -4A-16. 



(6.7) 



These give the power of s required for invariance of the solutions, namely s~ 2 and s~ p respectively. 

From the perspective of the LPA and scalar field theory, these power law solutions are surprising 
because they are not the ones expected by dimensional reasoning. Given that the mass-dimension 
of the associated coupling g is 2A, this reasoning would lead us to expect that vq ~ R 2 ~ x . However 
this reasoning is justified when the quantum corrections can be neglected in the large field regime. 
Here in contrast we see that the asymptotic behaviour of eigen-operators in this regime is controlled 



by the physical scalar quantum corrections (2.16). 

The fact that one of the solutions, namely v$ = R 2 , is independent of A, is also novel. The 



reason for this is that it solves separately both the term in curly brackets in (6.6) and the rest of 
this equation, i.e. to leading order separately both the left hand side and the right hand side of 



(6.3) vanish. If this were to persist at sub-leading orders it would mean that there is a solution 



to (6.1) whose i-dependence is undetermined since Sf = g(t)vo would solve the left hand side 
and separately 5f = g(t)vo would solve the right hand side, for any choice of g(t). We will see 
however that at sub- leading orders the two sides are re-coupled and furthermore, all three linearly 
independent solutions are coupled together by the constraints stemming from the finite singular 
points. 
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Expanding to next order, we look for a v\(R) that is ~ 0(1/ R) down (the precise definition 
being phrased as in (5.7) for Rvi). The next order comes from O(e ) pieces which involve again 



the third terms on each side of (6.3) and now also the first terms, originally from the left hand 



side of (3.2). The first terms just receive the now known vq- In the third terms, d(R) appears as a 
factor which cancels on top and bottom just as in the derivation of f\. The result is 



R 3 v'" - R 2 v" + 6Rv[ - 10vi + A {R 2 v'{ - Rv[ } = R^ 1 (a + c 2 cos In R 2 + c 3 sin In R 2 



(6.8) 



for vq = R p , and of the same form for vq = R 2 but with p = 2 on the right hand side. The c« are 



functions of A and the parameters A, B, C . Comparing (6.6), we see that the general solution just 
reproduces the three vq solutions and is therefore rejected. It is as before the special solution we 



are after. For the two powers (6.7) it is simplest to express the result in terms of p, inverting (6.7) 
to express A as: 

A = 2 — p — 5/p . 



(6.9) 



Combined with the leading order and setting e = 1 we thus obtain: 
_ x r3 (384ttV - 2p - 2h}A + p[p + 30]) 



v(R) = R P + 



R p - X i 



+ 



(P~ 3)(p 2 -p-5) 
768vr 2 (p 2 - 2p- 5) [t(-C, B) coslni? 2 + ^(B, C) sinlni? 2 ] 



(p 2 -6p+ 13) (p 4 - 2p 3 - hp 2 + 10p + 25) 



} + 0(R p - 2 ) 



where we have written 



7(2, y) := 2(2p 2 - 4p - 5)x + (p 3 - 4p 2 - 6p + 15)y . 



And for the remaining solution we obtain: 

'6(960vr 2 Ayl - 6A + 1) 768vr 2 



v(R) = R Z + R 



A + 4 



+ 



([B - 2C] cos In R 2 + [C + 2B\ sin In R 2 



(6.10) 



(6.11) 



+ 0(1). 
(6.12) 



Again this is as far as we can take the asymptotic expansion in closed form. The next and higher 
orders can be expressed as an integral of ratios of terms, similar to the asymptotic expansion for 



f{R)- Since (6.3) is linear in v, the denominator terms are only those generated by /. Therefore 



we obtain no further asymptotic restriction on parameters than (5.17). 



So far we have assumed that A is some generic value. There are a host of special cases. In (6.6), 
if A = — 5/2 or A = 2 ± \/5 then the powers are degenerate, since in the first case p = 2, |, and in 
the second two cases the square- root vanishes in (6.7) and the two values of p coincide as p = 
respectively. Thus in these cases there is a further solution R p lnR 2 , where p = 2 or p = 
respectively. Higher order terms for these new solutions can be developed as above. 



30 



If A = —4 in (6.12), or such that it corresponds to p being a root of one of the denominators 



in (6.10), then these expansions are no longer applicable. Instead the special solutions need to be 
recomputed and have a new piece proportional to RlnR 2 . After this the higher order terms again 
may be computed as above. 

We will ignore these special cases in what follows. The main property we will require is that 
there are three unique asymptotic solutions, i.e. apart from an overall scale factor, fixed here so 



that the first coefficient in (6.10) and (6.12) is set to one, there are no further parameters in them 



beyond A and the A,B,C parameters inherited from f(R). This is true for the special cases also. 



Using these as asymptotic boundary conditions to (6.3), we thus find three unique solutions, which 



however can be expected to develop singularities at the fixed singular points. 

Since polynomial truncations result in complex eigenvalues |2T[6| , and there is no obvious reason 



for (6.3) to exclude complex eigenvalues, we will not assume that A has to be real. In (6.3), the 



coefficients of v ^ and Xv ^ are real. Therefore, if a complex A and thus complex v, satisfies the 
equation, then so does A and v. In this case we obtain the required real perturbation by adding 



the two t-dependent complex conjugate eigen-perturbations (6.2) together with complex conjugate 
coefficients a. 



If A is complex, then from (6.9), evidently p is also complex. Since (6.9) is invariant under 



p \— > 5/p, if p is one power corresponding to A, then the other power is 5/p. From here on, we use 
this freedom to define p such that \p\ > y/5. If \p\ > v5 then this picks out uniquely one of the 
roots. If |p| = then the other root is p. In this case, writing p = \/5expi#, we see that 



A = 2 - 2\/5cos6> 



(6.13) 



is real. In this case the two asymptotic solutions (6.10) should be added together with complex 



conjugate coefficients, yielding a real result where the complex powers can be re-expressed in terms 
of i^ cose cos(V5sin01ni?) and R^ cosd sm(V5sm9lnR). 

Define now w+ to be the unique solution whose behaviour as R — > oo satisfies (6.10) for our 



chosen p. Similarly define W- to be the unique solution whose behaviour matches (6.10) for power 



p i — y 5/p. Finally define w to be the unique solution that matches (6.12) as R — > oo. The general 
solution for v (R) is therefore given by: 

v(R) = a + w + (R) +a-W-(R) + w(R) (6.14) 

where the a± are two coefficients, and for the sake of argument we have chosen to normalise v(R) so 



that as R — > oo, its R term has unit coefficient, as follows from (6.12). (This is allowed providing 
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this coefficient does not vanish. The following arguments can be generalised as appropriate if this 
happens.) 

For later use we note that if Re(p) > 2, w+(R) will be the dominant contribution at large R, 
so that v(R) rj a+RP in this regime, and from ( |6.9[ ), 

ReA = 2 - Re(p) - 5Re(p)/\p\ 2 < -10/|p| 2 , (6.15) 

hence these correspond to irrelevant operators. 

On the other hand if Re(p) < 2, then w(R) is the dominant contribution at large R, implying 
that in this regime v(R) ~ R 2 . In this case, ReA > — 10/|p| 2 , so the relevancy or irrelevancy of 
these operators is not yet determined. 

6.2 Counting parameters in the eigen-operator solutions 

If A is real and satisfies 

|A — 2| > 2V5 , (6.16) 
then p is real and the a± give us two real parameters. If A is real but does not satisfy the 



above, then it is of form (6.13) and the a± are a complex conjugate pair. Therefore again the 
a± supply two real parameters. If A is complex, then the a± supply two complex parameters. In 
summary, in the appropriate field over which A is defined (viz. real or complex), the ct± supply two 
parameters. Therefore so far we have no restriction on the number of parameters in the solution 
v(R) (or their domain) beyond the fact that we can eliminate one of them with a normalisation 
condition. As we have already mentioned in the introduction, in the LPA for scalar field theory the 
asymptotic behaviour of the eigen-operators divides into two classes: those that grow as a power of 
the field ip and those that grow as the exponential of a power. The latter are excluded by several 



arguments 22 23 , as already reviewed at the end of the introduction. Here there seems to be no 



reason to exclude some of the asymptotic behaviour in (6.10) or (6.12). We will explore this in 
more detail in sec. 16.31 



As we found in the fixed point equation, there are no fixed singularities coming from the 



denominators in (6.3); they only arise from the coefficient of v . Since this is the same as the 



coefficient of / in (3.2), we have again three fixed singularities with locations R = R c = 0,R±. 
Since the differential equation is linear in v, requiring that v is smooth means that these single 
poles provide linear constraints on the derivatives of v at R c ; they take the form 

h 2 v"{R c ) + hi v'(R c ) + ho v(R c ) = . (6.17) 
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As is obvious from (6.3), the coefficients have a linear dependence on A as hi = h]\ + h®, and the h\ 
have complicated dependence on f(R c ),f'(R c ) and f"(R c ). Unfortunately even the constraint at 
R = is too long to make it worthwhile displaying. The main point is that as a result, a normalised 



smooth solution v(R) satisfies a constraint (6.17) at each singular point. 

If we take the view that the applicable domain is < R < oo, the first of the situations discussed 
in sec. 



4.3 then only the singular points at R = and R = R+ are operative. For the same reasons 



as discussed in sec. 4.2 we can be certain that the corresponding constraints act independently. 



Therefore for each A, the two parameters a± in (6.14) are fixed to a discrete set of values. Regarding 
A as just another parameter, we have three parameters but only two constraints, therefore we will 
find solutions for continuous range(s) of A (two-dimensional domain(s) if A is complex). Now, given 
that the fixed point solution f(R) is one of a line of fixed points, we know that one eigen-value is 
A = 0, corresponding to the perturbation that takes us to a new fixed point f(R) + Sf(R). There 
is no reason to expect this eigen-value to be always at the end of a range, therefore typically the 
neighbourhood of A = also provides solutions, and thus includes Re(A) > 0. Thus we see that 
almost always, each fixed point on the line supports an uncountably infinite number of relevant 
eigenvalues. Since these would correspond to an uncountably infinite number of renormalised 
couplings, this is of course not an acceptable starting point for a predictive theory of gravity. The 
only way to escape this conclusion is if the A = eigen-value lies at the end of a range. In this 
case however we will either again find that Re(A) > provides an uncountably infinite number of 
relevant eigen-values, or Re(A) > is excluded - in which case we have no relevant operators at 
all and thus in the limit the theory has no adjustable couplings and is scale- invariant. As already 



discussed in the introduction and sec. |4.3[ the problem is that the eigen-operator equation is not 
yet over-constrained. 



As we discussed in the introduction and sec. 4.3, this situation is resolved if we smoothly match 



into spaces with R < 0, as we argue we physically are required to do in any case. Taking (3.2) 
as working example, we now have all three singular points (R = 0,R±) and these supply three 
independent constraints. Fixed points themselves now form a discrete set. As discussed in sec. 
[7J unfortunately for the current equation (3.2) this set is almost certainly empty. If there were 
such a fixed point then, given A, the two parameters a± are over constrained. Regarding A as 
one more parameter, we see that we would now find a discrete set of solutions and thus quantised 
renormalisation group eigen-values, as expected. 
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6.3 Asymptotic properties of finite perturbations 



As already mentioned, cf. below (6.7), the powers involved in the asymptotic solutions (6.10) and 



(6.12) are not the ones expected from dimensional reasoning. This can be traced to the fact that 
the quantum corrections in the flow equation do not decouple for large R (see also the comments 
in the conclusions). On the other hand, in the analysis of the LPA for scalar field theory at large 
field, the recovery of the powers expected by dimensional reasoning is an important first step in 
proving that the eigen-perturbations have sensible properties when used as small but finite seeds 
for the renormalisation group flow out of the fixed point, distinguishing these from the non-power- 
law perturbations which are superfluous and should be discarded, as already reviewed at the end 



of the introduction, also see [22,23 . The question then arises whether all of the solutions (6.14) 
should be kept when viewed in this way or whether some of these should also be discarded. We 
have already seen that if R < is consistently included, the eigen-perturbations form the expected 
discrete spectrum, and we will see that the asymptotic regime paints a self-consistent picture for 
all of these operators when treated as finite perturbations. However the analytic study of the large 



field regime of the present equation (6.1 ) is not as powerful as it is in the case of scalar field theory, 
the ^-evolution ultimately being controlled by the finite singular points. 



The separable solution (6.2) is valid providing a is small enough to ensure 

5f(R,t) 



(6.18) 



so that the linearisation step in (6.1) is valid. 



We first consider the status of (marginally) irrelevant operators with small but finite couplings. 
These thus have ReA < 0. At the end of sec. |6.1| we saw that their asymptotic behaviour depends 
on whether ReA is greater or less than — 10/|p| 2 . 



If we assume that 13 ReA > — 10/|p| 2 , then recalling (4.14), both the perturbation and the fixed 
point solution grow like R 2 at large R. Therefore, since these solutions are smooth, if (6.18) is valid 
for some finite R it will be valid for all R. Since the t dependence is given by a shrinking exp — 2\t 



(or a more general shrinking function if the operator is marginally irrelevant) if (6.18) is valid at 
the overall cutoff scale t cu t = In k cu t/n then it is valid for all t as it is lowered below the cutoff scale. 
We conclude that these operators correspond to genuine irrelevant perturbations that fall back into 
the fixed point following the already determined t-dependence, as the cutoff k is lowered. 

Now consider those cases where ReA < — 10/|p| 2 . The perturbation v(R) now grows faster, as 
RP where Re(p) > 2. Therefore for any finite choice of a and t, there is some R beyond which 



If A is real one can show that this corresponds to A > — 2. 
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the separable solution (6.2) no longer satisfies (6.18). If we substitute f(R,t) = f(R) + 5f(R,t 



directly into (2.12) and use the fact that 5f(R,t) ~ R p dominates at large R, we see that actually 
in this case the quantum corrections on the right hand side can be neglected (growing no faster 



than R ) and thus 5f(R,t) has to satisfy the left hand side of (2.12) on its own. Therefore in this 



regime, 5f(R,t) again takes the separable form (6.2) but with A now taking the value expected on 



dimensional grounds, namely A = 2 — p. Comparing with (6.9), we see that as we lower t, the large 
R regime is controlled by an effective renormalisation group eigen- value A = 2 — p which is still 
irrelevant, but less so than the one derived at small R and for truly infinitesimal perturbations. In 
summary, although some irrelevant perturbations, when considered small but finite, behave in an 
unusual way, all of them still fall into the fixed point as we integrate out modes by lowering t. 

We can therefore ignore the irrelevant perturbations now and consider what happens when we 
build the renormalised trajectory by perturbing by some small but finite linear combination of the 
(marginally) relevant operators. In this case each v(R) ~ R 2 for large R, therefore together we find 



that in this domain, Sf(R,t) = g(t)R for some finite function g(t). If (6.18) is valid at finite R it 



will thus also be valid as R — > oo. As we noted below (6.7), at the linearised level g(t) is actually 



not determined by constraints in this regime. Similarly to the fixed point, cf. the discussion around 



(5.12), even for the full equations (2.12) this limit is singular and makes sense only when developed 



to higher order. At the linearised level g(t) is determined by the matching at the singular points 



R = 0, R± which couples together the other asymptotic solutions in (6.14), as explained in sec. 6.2 
Thus we learn that actually g(t) must be expressed as 

g (t) = ^ a„ exp -2A n t , (6.19) 

n 

where the sum is over the (marginally) relevant operators. As modes are integrated out, t falls from 



t cu t, and g(t) grows. For some t <C t cu t, (6.18) will no longer be satisfied, however this will happen 
at the singular points and elsewhere at more or less the same t as it happens in the asymptotic 



regime. From here on the t dependence changes from (6.19) as determined by the full non-linear 



equation (2.12) and in particular by a full non-linear matching at the finite singular points. 



7 Numerical solutions 



7.1 General approach 



The numerical analysis of the fixed point equation (3.2) is dominated by the presence of the three 



fixed singularities i?+,0,i?_. It is therefore sensible to divide the real axis into the following 
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intervals 

J_ 00 = (-oo,i?_], J_ = [fl_,0], l+ = [0,i?+], /oo = [i?+,oo). (7.1) 

A careful numerical search for solutions has to be carried out separately on each interval and, as 
mentioned before, any possible solutions for R < have to be regarded as tentative solutions since 



(3.2) was derived under the assumption R > 0. In order to bridge the fixed singularities we Taylor 
expand any potential solution around them. In the case of R+ we have 

f(R) = b + b 1 (R- R + ) + V b -^hhl(R - R +T . (7.2) 

n=2 

Similarly, 

f(R) = a + a 1 R + T °" (q °' Ql } R" (7.3) 

Z — / n \ 



n=2 



is the Taylor expansion around and finally 

5 



^— ' n! 

n=2 



/(i?) = Aj + /3i(JZ - + V Pwl - ^-) n - (7-4) 



bridges across 

In each case, the requirement that the solutions f(R) be regular at the singularity (cf. sec. [4]) 
translates into the fact that only the first two coefficients are independent. All higher coefficients, 
starting with the second, can then be expressed as functions of the first two. It should be noted 
that these expressions are unique for all coefficients except 62 and 02 for which we find quadratic 
constraints. Their solutions can be uniquely written as 

bt(bo,h) = (p(h,h) ± (7-5) 

and 

(3f(f3 ,(3i) = -726ft- (p(Ai>Pi) ± yMfafc)) > (7-6) 

where p, q and p, q are polynomials in bo, b\ and /?o, (3i, respectively. As we will see, both possibil- 
ities for 62 and (5% lead to fixed point solutions. As displayed above, the Taylor expansions we have 
used for our computations are of order 5 with the higher coefficients quickly becoming complicated 
expressions of the first twoj^] Especially for the expansions around R± it would be a challenging 
task to handle them for even higher coefficients. 



In fact we extend (7.31 to sixth order for part of the analysis for reasons discussed at the end of sec. 



7.4 
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Our search for numerical solutions was performed within the Maple(TM) package |34| employing 
the shooting method. We picked the numerical integration method "dverk78" which is a seventh- 
eighth order Runge-Kutta integrator and is able to operate at arbitrarily high precision. 

We now illustrate the individual steps that led to the discovery of solutions starting with the 
region R > 0. If we started shooting down from some large R&Q to R+ we would face the problem 



of having to search through the 3-dimensional parameter space of the asymptotic expansion (4.14). 
It is therefore more reasonable to shoot out from R + or where the parameter space is only 
two-dimensional. Moreover, we start looking for solutions on I + as opposed to shooting from R + 
towards large R as we expect the additional constraint at to reduce our parameter space by one 



dimension which is a more severe condition than (5.17) for the asymptotic series. Most solutions 
we were able to find can be obtained by both shooting out from zero towards i?+ and from R + 
towards zero, in particular those which are valid for all R > 0. The exceptions were solutions valid 
only on I + with a moveable singularity appearing shortly after R + , e.g. at R c ~ 8; in this case 



integrating from R + to zero turned out to be difficult 5 For the purposes of illustration however 



let us focus on shooting out from zero towards R+ in the following. 



Given any initial pair (ao,ai) we use (7.3) to compute the three initial values 



f(eo), f'(e ), f"(e ) (7.7) 

at distance £q to the right of zero. Shooting out from zero towards R+ we will find that most pairs 
(cjq, oi) end at a moveable singularity. For some however, the numerical integrator supplies us with 



the solution f n valid close enough to R + so that by using (7.2) we can solve the system 



f(R+ - ei) = f n (R+ - £i), f'(R+ - ex) = f n (R + - ex) (7.8) 



at B\ away from R + for bo and b\. Using (7.2) again, we are then in a position to compute 



f"(R+ — £i) and compare it to the actual value f%(R+ — £x) via the quantity 

5f ^ = nR + -ex) • (7 - 9) 

This represents the additional matching condition at R+ and implements the constraint coming 
from the singularity at R = R+ . A solution is found if the two second derivatives agree to sufficient 
accuracy. In fact, our criterion for a valid solution on /+ will be that 5f so \ varies smoothly across 



zero upon variation of ax while keeping ao fixed. This condition on 5f so \ is illustrated in fig. 7.1 



Small variations of ax lead to the expected linear behaviour of 5f so \. Varying ax by a larger amount 



3 Similarly, integrating from R- to zero was ill behaved for partial solutions with a moveable singularity at R c 
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Figure 7.1: The quantity Sf ao \ varies linearly across zero as a function of a\ keeping ao fixed. If 
this is the case we recognise the pair (ao, ai) as a solution point for I + . 



brings the non-linear form of Sf so \(ai) to light whereas if we go to very small variations of a\ we 
find numerical fluctuations superimposed on the linear behaviour. There are good reasons for this 
solution criterion and why we should not just require |5/ so i| to be small, as we will explain in detail 
in sec. 17.41 

A priori there is no special region in the a-plane of points (ao,ai) where solutions can be 
found. However, there is an important structure to the a-plane which is connected to the Taylor 



expansion (7.3). The coefficients 02(00, ai), 0,3(0,0, a\), . . . become singular along the dashed lines 



7i in fig. 7.2 which implies that the expansion (7.3) is not valid anymore. For reasons that will be 



explained in sec. 7.4 along these lines there are no solutions except at finitely many points. It is 
therefore reasonable to search for the first solution between these singular lines. We can do this by 



using similar graphs as in fig. 7.1 where we plot 5f so \(ai) m the range between two singular lines 
while keeping ao fixed or, if moveable singularities prevent the solver from reaching _R+, instead 
of o/soi( a i) w e can plot the maximum R reached in a similar manner. In this way it is possible 
to find the first solution point in the a-plane. Once we have a solution (ao,ai), we create two 
new a-vectors (ao + oao, a\ + 5a\) and (ao + 5oq, a\ + <5ai), where the two variations (5oq, 5a\) and 
(oao, 5a\) are orthogonal to each other and small enough so that we are still able to integrate from 
zero to R + . Generically 5f so \ will be large at these new points and so we repeatedly apply the 



secant method to 5f so \(ao, a\) until we find a new pair (ao, ai) with the behaviour shown in fig. 7.1 
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We have automated this process by implementing it in a procedure in Maple and it is this method 
of obtaining a new solution point in the a-plane starting from an existing one that lies at the heart 
of how we found our solution lines. Given at least two solution points (ao,oi) and (So,5i) we can 
optimise the procedure we have just described by using their difference as a guide for the direction 
in which to look for the next solution point, i.e. we define 

(Sao, Sax) oc (clq — ao, ai — ai). (7.10) 

These two methods constitute the way by which we extended into lines of solutions starting with 
only one solution point. 

Clearly, the numerical approach to the search for solutions on J_ proceeds in a similar way. 



We start off with a pair (ao,ai) and use (7.3) at eo to the left of zero, to determine the initial 



conditions which are needed for shooting towards R-. If the (partial) solution makes it to R- + £1, 



we solve the corresponding version of (7.8) here and adapt ( |7.9| ) to judge whether we are dealing 
with a solution. We can employ the same methods as for /+ to find the very first solution on J_ 
and to extend into solution lines. 

When it comes to finding solutions on Iqq the degree of difficulty of the problem at hand has 
already been lowered by the fact that we can restrict ourselves to check if any of the solutions 
valid on i+ will make it out to infinity, i.e. we are left with a one dimensional parameter space. 



The strategy therefore is to select a possible candidate pair (bo, b\) and to exploit (7.2) at i? + + E\ 



to find the initial values. We then integrate up to some large Roo where we finally match to the 



asymptotic expansion (4.14). The accuracy of the three asymptotic parameters A, B, C will depend 



on the value of R^ which is why it should be taken reasonably large] 16 If the constraint (15.17) 



is satisfied the existence of the solution is guaranteed for all R > Roq. On the other hand, if the 



asymptotic parameters A,B,C are such that (5.17) is not fulfilled, we are dealing with a partial 



solution only since a moveable singularity is bound to appear at some R > Roo, cf. sec. [5j In 
principle this same method can be used to find solutions on I-oo but as we will see, none of the 
solutions we found, on 7_ extend, to U I—qq. 

7.2 The solutions 



After extensive numerical searches using the approach outlined in sec. 7.1 we were able to obtain 
a total of five solution lines valid on the range R > 0, four additional solution lines representing 



3 We defer a discussion of this important point to sec. 7.3 
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solutions on J_, and we found evidence for more solution lines. For clarity of presentation we devote 
a separate section to each set of solution lines. 



7.2.1 Solutions for R > 

The lines of solutions defined on the original domain of validity I + U 1^ of the fixed point equation 



(3.2) are shown in fig. 7.2 We have three lines of fixed point solutions in the range ao > 0, all 
starting at the origin and moving out between two singular lines of the expansion ( |7.3[ ) f^j In fact, 
the bottom two solution lines are both situated between the two singular lines 73 and 74 in the 



region close to the origin, cf. fig. 7.5, but the black solution line crosses the singular line 73 in its 



early stages and moves out between 72 and 73. In order to find these fixed point solutions in the 



region of positive ao it was crucial to use the 63 (&o>&l) solution in (7.5). Furthermore, all these 
solutions are valid for all R > 0, i.e. satisfy the asymptotic constraint (5.17), and their coefficient A 
of the dominant asymptotic component is positive. As explained in the Discussion and Conclusions, 
A > is a necessary (but not sufficient) condition for the solutions to make physical sense. The 



corresponding solution lines in the o-plane are shown in fig. 7.3 (la), (lb). All three start at the 
origin and approximately follow a straight line towards larger values of 60 • Even after subtraction 
of this linear component, the result being shown in fig. |7.3|(lb), it is impossible to distinguish them 



by eye. The distance between these solution lines varies but can be as small as Ab\/b\ = 10 



-10 



In the range ao < there are two more solution lines but they behave differently, see fig. 7.2 
Leaving the origin they run towards the left but eventually both turn around to approach the 
origin again. In contrast to the solutions with ao > 0, these solutions all have A < 0, i.e. they 
assume negative values for large enough R. In order to satisfy the constraint arising from the fixed 



singularity at R + it was necessary to employ the 6^"(6o 5 0i) solution in (7.5) in this region of the 
a-plane. If we look at the corresponding solution lines in the 6-plane shown in fig. 7.3 (2a), (2b) 
we again observe a linear behaviour. Unlike in the case for 60 > however, we can distinguish the 
two solution lines clearly after subtracting their straight line approximation. It then also becomes 
clear that, just like the solution lines in the a-plane, after moving away from the origin they turn 
around and approach it again. In both the a- and the 6-plane the blue solution line is contained 
inside the region bounded by the red solution line. 

The fact that the solution lines in the 6-plane can be approximately described by a straight line 



Further such solution lines are expected between singular lines j n with n > 5, as we discuss in sec. 7.4 
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Figure 7.2: Five lines of solutions in the a-plane corresponding to the Mines in fig. 7.3 Different 



solution lines are represented consistently in different colours (these are red and blue in the range 
ao < and, from bottom to top, black, magenta and green in the range ao > 0). The black points 



represent example solutions which are plotted in fig. 7.4 Along each of the black dashed lines a 



coefficient in (7.3) becomes singular (c/. sec. 7.4). 



can be attributed to the following observation. If we rescale according to 
f(R) = bJ(R) =bi(^ + (R- R+) + b f i (R~ R+f + 



and consider the limit b\ — > 0, the fixed-point equation (3.2) becomes 



r 2 + 7i + v p + % 



f=f 



(7.11) 



(7.12) 



This follows from the fact that the right hand side in (3.2) is invariant under rescaling of / but 



the left hand side is linear. Hence we would expect that, close to the origin of the 6-plane, any 
solution line should show predominantly linear behaviour. The same conclusions can be drawn for 
solution lines in the a-plane. However, computing the slope of the solution lines in this plane as we 
approach the origin, we find the slopes vary linearly, implying quadratic behaviour of the solution 
lines close to the origin. Given the huge magnification factor implied above, which occurs as we 
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Figure 7.3: Three lines of solutions on 7+ U loo in the range bo > in (la) and two additional 
solution lines in the range bo < in (2a). In both cases we have subtracted the straight line given 
by b± = 0.254352 • bo — 2.5 • 1CT 7 and plotted the result on the right. 



integrate from R + to zero, we would expect to recover linear behaviour of the solution lines at 
distances much closer to the origin of the a-plane than we have probed here. 



The non-linear components of the solution lines in the 6-plane in fig. 7.3 where obtained by 



subtracting the same straight line from both the solution lines with positive and negative bo- This 
is a first indication that it could be possible to match these solution lines across the origin. If one 



computes the second coefficient 62 along the solution lines the two possibilites (7.5) also seem to 



match across the origin. However, since the Mines are so close to each other it is difficult to give a 
definite answer to this question in the 6-plane. The solution lines in the a-plane can be separated 
clearly and computing their slopes close to the origin it seems to be possible to match the lines 
with positive ao uniquely to a solution line with negative ao, e.g. the black solution line matches 
to the upper part of the red solution line and the solution line in magenta can be matched to the 
upper part of the blue solution line. The lower part of the blue solution line can be continued 
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by the green solution line whereas the lower part of the red solution line seems to belong to an 
as yet undeveloped line on the other side of the origin. Although these pairings are confirmed by 



matching the other coefficients of the a-series (7.3) across the origin, none of these methods would 



be able to distinguish solution lines that become tangential as they approach the origin. 



The solution lines on the right hand side in fig. 7.2 can be continued towards larger values of 
do but the numerical integration becomes ever more time consuming. We have also found a first 
solution point on four additional solution lines between 74 and 75, on both sides of the origin. 
Judging by these first solutions they consist of partial solutions only. They could also complicate 
the matching of solution lines across the origin since three of them are very close to each other 
compared to the solution lines we know of and could thus be tangential at the origin. 



As indicated by the points in fig. 7.2, we have selected one example solution on each of the 
solution lines for R > with A > and plotted the fixed point function f(R) and its first and 
second derivatives in fig. |7.4| All three fixed point functions assume positive values only. The 
first and second derivatives already show that we approach quadratic behaviour close to i?+ as the 
plots on the right hand side in fig. |7.4| confirm. Nevertheless, in agreement with what is already 



built into the asymptotic expansion (4.14), we can also discern the oscillatory pattern from the 
plots of the second derivative. The parameter values for these solutions are listed in table [TJ In the 
Discussion and Conclusions, we comment on the physical interpretation of these solutions and also 
the solutions with do < 0. 



FP 


a [10- 4 ] 


ax [10- 4 ] 


bo [10- 3 ] 


61 [1CT 3 ] 


A [10" 5 ] 


B [10" 5 ] 


C [10" 5 ] 


green 


6.1053 


8.6648 


21.869 


4.8337 


28.741 


-3.2321 


-1.3074 


magenta 


5.4421 


5.7113 


14.899 


3.2125 


19.207 


-3.3397 


-0.47817 


black 


8.5841 


8.6461 


42.986 


10.002 


59.479 


-1.5161 


-2.1547 



Table 1: Parameter values for the fixed point functions plotted in fig. 7.4 Note that many of these 



parameters are known to higher accuracy than given here, cf. sec. 7.3 and table [2l 



7.2.2 Solutions for < R < 



The fact that the fixed point equation (3.2) admits the uncountably infinite number of solutions 
presented in the previous section can be expected from the parameter counting method. Although 



the derivation of (3.2) was carried out with spaces of positive scalar curvature R we have neverthe- 
less extended our analysis to the range R < in the hope this would constrain the solution space 
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Figure 7.4: Fixed point functions and their first and second derivatives represented by the points 



shown in fig. 7.2 in corresponding colours (from bottom to top these are magenta, green, black) 



plotted over the interval I + in the left column and over the range R + < R < Roo on the right. 

to at most a finite number of physically sensible fixed point functions valid for — oo < R < oo. 
Since in this regard we are interested in fixed point solutions with positive quadratic asymptotic 
component, i.e. A > 0, we have concentrated our efforts on finding solutions valid on / which are 
situated in the ao > half-plane. 

We have been able to find four solution lines in this range whose form in the a-plane is displayed 



in fig. 7.5 Two of these solution lines extend as far as the previous solution lines for R > whereas 
the other two turn around while still close to the origin running alongside the bottom singular line 
72. Because the latter two lines are so close to 72 we have refrained from developing them all the 
way back to the origin but we have checked that they are still present at ao = 10 -4 where they are 
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approximately 2 • 10 9 away from 72. 



The situation is different for the brown and olive solution line in fig. 7.5 They start out between 
the same singular lines 73 and 74 as the two solution lines for positive R shown in black and magenta 
and both cross 73 where the black solution line does but approach the black solution line for large 
values of ao- The olive solution line then intersects the black solution line and thus gives rise to a 
solution valid on the interval R < R < 00. Although the brown solution line does not cross the 
black solution line again in the range we have developed it in, it is to be expected that it will do 
so if the two lines are continued further. Close to the origin we find two more solutions valid on 
/?_ < R < 00: the solution lines in brown and olive both intersect the solution line in magenta for 
I + exactly once (this is barely discernible in fig. 7.5). The former at (ao, a±) = (9.8- 10 -5 , 8.1 • 10~ 5 ) 
and the latter at (ao, ax) = (6.07 • 10~ 5 , 4.9 • 10" 5 ). 
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Figure 7.5: Four solution lines for J_ in the a-plane together with three of the solution lines 



displayed in fig. 7.2 on the left and a magnification of the region close to the origin on the right. 



In the magnification the colouring from bottom to top for the solution lines on 7_ is orange, violet, 



olive and brown. The black dashed lines are the same singular lines 7, as in fig. 7.2 



The corresponding solution lines in the /3-plane are shown in fig. 7.6 Like the solution lines 
in the fo-plane they are governed by an underlying linear behaviour with superimposed non-linear 
variations which in the case of the left plot are of order 10 -8 and for the plot on the right they 
can be as large as 3 • 10~ 5 . Note the two different ways in which these four lines are paired: in the 
a-plane the two long and the two short solution lines are paired together whereas in the /3-plane 
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it is always one /3-line corresponding to a long solution line in the o-plane that is paired with a 
/3-line corresponding to a short solution line in the a-plane. As we would expect from continuity, 
varying the initial pair (ao,oi) along a smooth curve from one of the short solution lines in the 
a-plane to the other leads to a smooth curve of points (/3o,/3i) in the /3-plane connecting the two 
corresponding solution lines in the two different regions of the /3-plane. However, this connecting 
curve does not consist of solutions since the constraint at R- is not satisfied. Curiously, despite 



0. 1093- 



0.1092- 



0.1091 



0.1090- 



0.1089- 



0.1088- 




0.0770- 



0.0768- 



0.0766- 



0.0764- 



0.0762- 



0.0760- 



-0.1158 -0.1154 ' -0.1150 



0.1144 -0.1140 




-0.075 -0.074 -0.073 -0.072 -0.071 -0.070 



Figure 7.6: The four solution lines for I_ in the /3-plane in the same colours as the corresponding 



lines shown in fig. 7.5 



the fact that the solution lines in the a-plane approach the origin, the corresponding lines in the 
/3-plane do not. Another difference to the situation for R > is that we have not found solution 
lines between the upper two singular lines 74 and 75. 

We have also found evidence for a line of solutions valid on J_ between the middle two singular 
lines 73 and 74 in the range of negative an- As we would have to match any such solution at R = 
to a solution with negative asymptotic parameter A, these solutions are physically less relevant. 

It has to be emphasised that none of the solutions for /_ presented here extend to /_oo U J_ . 
Instead they all end at a moveable singularity R c ~ 8. 

7.3 Error analysis 

The actual numerical computations have been carried out with 25 significant digits in Maple and 



we used the values e$ = 5 • 10 and £\ = 10 in (7.7) and (7.8) and also in the matching for 
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R < 0. We will now go through a detailed error analysis for the example solutions plotted in fig. 



17^1 



The first source of error in our solution strategy comes from the truncated Taylor expansions 



(7.2)-(7.4). We investigated the radius of convergence of (7.3) and (7.2) by performing a root test 



on the first 60 coefficients for each of the example solutions. The result for the green example 



solution is shown in fig. 7.7 The root test on the left hand side points to a convergence radius 



Ml " 3 




Figure 7.7: The result of a root test applied to the first 60 coefficients of (7.3) on the left and of 



(7.2) on the right for the green fixed point solution in fig. 7.4 



of ro ~ 0.4 for the Taylor series around zero. The plot on the right pertains to a root test for the 
Taylor series around i?+ and seems to indicate that the convergence radius of that series is r + ~ 9. 
The two convergence radii have to satisfy the relation r + < R + + tq. Since R+ ~ 7.4, this is roughly 
satisfied for the above values for ro and r + . 

If we use these Taylor approximations and plot them together with the numerically computed 
solution we get fig. |7.8| We notice the rapid divergence of the Taylor expansions from the numerical 
solution towards the end of the respective convergence interval but in between the numerical solution 
is indistinguishable from the Taylor approximation. 



Zooming in on the Taylor approximations in the region around R = 0, cf. fig. 7.8 (right), we 
can see that for negative values of R the 6-series approximation in blue starts diverging earlier than 
the red o-series approximation. This plot is consistent with a radius of convergence r + ~ R + + 0.4 
at R = i?+ which we would need to satisfy the aforementioned relation between r + and rg. 



The analogous plots to fig. 7.8 for the first and second derivatives of f{R) for the green solution 
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Figure 7.8: The green example solution together with its Taylor approximations of order 60 around 
i?+ and zero in blue and red, respectively. The magnification on the right shows their behaviour 
round R = 0. 

example show equally good agreement with the only difference that the Taylor expansions start 
diverging slightly earlier as the number of derivatives increases. Altogether these results present 
strong evidence that we can trust the Taylor expansions used to bridge the singularities. 

With this in mind, the error analysis proceeds along the following lines. We estimate the 



relative error A ajS in truncating the Taylor series (7.3), representing the relative error on the initial 



conditions (7.7), by using the modulus of the last term in the Taylor expansion. If we now suppose 



the true initial conditions (7.7) are within this error of the ones we actually have, we can determine 
the amount by which we have to vary the two parameters (do, oi) in order to obtain this variation 
A a s of the initial conditions and then take the result A Q)S as the corresponding error on the pair 
(a ,ai). 

We also have to consider the intrinsic error of the solver, i.e. even if we started with exact initial 
conditions at Eq, the integrator will induce a small error in each step which leads to a numerical 
solution which is different from the true solution. The size of this error A a ^ n can be determined by 



using plots like fig. 7.1 and varying a\ by ever smaller amounts. At some point we will be able to see 
numerical fluctuations appearing and we have reached the maximum accuracy. These fluctuations 
limit the accuracy to which we can determine the zero of <5/ so i we are looking for. The value of 
n can then be read off from the ai-axis of this plot. Since these two errors are uncorrelated we 
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obtain the total error on (ao,a±) according to 



A a = JAL + A2 . (7.13) 



To obtain the resulting error on the corresponding solution point (60,61) we now evolve this total 
error A a with the numerical integrator by starting from the modified pair (1 + A a ) • (ao,ai) and 
compute the error A^ in on the pair (6q, 61). We then also have to include the truncation error of the 



Taylor series at R + , viz. (7.2), in a similar manner as for the truncation error at zero, and finally 



compute the total error A^ on the pair (60,61) with the corresponding form of (7.13). 

By similar methods we can assess the error on the asymptotic parameters A, B, C. At R^ we 
have two sources of errors. Firstly, we have to take into account that the error on the point in the 
6-plane will have evolved to a new error Aoo jn on the initial conditions at R^ and hence on the 
asymptotic parameters A,B,C. Secondly, we also have a truncation error Aqq s produced by the 



asymptotic series (4.14). Again, we combine these two errors into the total error on any of 



the parameters A, B, C by using the analogous form of (7.13). Compared to the errors on solution 
points in the a- and 6-plane the error on the asymptotic parameters is relatively large. This is 
to be expected from the fact that we have to integrate over a comparatively long range, namely 
from i?+ to Roo = 200, in order to determine these parameters and thus the error on the solution 
point in the 6-plane can evolve to larger values. Additionally, we truncate the asymptotic series 
at order R°, i.e. we keep only the first three orders, which also contributes to the size of A^. 
Despite the magnitude of Aqo, the essential message delivered by the asymptotic parameters, that 



the constraint (5.17) is fulfilled, remains unaltered for all example solutions. 

The various values for the errors on points in the a- and 6-plane as well as on the asymptotic 
parameters are summarised in table [2} 

We have verified that the position of these example solutions in the a- and 6-plane remains 
unaltered if we decrease the distances Eq and E\ used to compute the initial conditions close to zero 
and This can be used to shrink the series errors A a>s and A^s in table [2] to smaller values. In 
fact, the value for Ab ]S in the case of the black example solution has been computed with e\ = 10 -4 
instead of the usual value E\ = 10 -3 which has been used in all other cases. In this way we can make 
sure that the different solution lines in the 6-plane are not an artefact of truncating the 6-series and 
would still be present if we used an arbitrarily long Taylor series at R + despite the fact that they 
can get very close to each other. 



49 



FP 


green 


magenta 


11 l 

black 




4 • 1CT 11 


3- 1(T 12 


3 • 1(T 14 




l • icr 14 


1 • 1CT 14 


1 • 1(T 12 


A 

A a 


4 • icr 11 


o i n— 1 2 

3-10 iZ 


i in— 12 

1-10 


A 6)S 


2.8 • 1CT 7 


1.4 • 1(T 7 


1 • 1CT 9 


A&, n 


3 • 1(T 9 


2 • 1CT 10 


3 • 1(T 10 


A 6 


2.8 • icr 7 


1.4 • 1(T 7 


1 • 1CT 9 


Aoo,s 


3.5 • 1(T 5 


5 • 1(T 5 


2 • 1(T 5 


Aoo,n 


1.8 • 1(T 5 


1.3 • 1(T 5 


3 • 1(T 7 


Aoo 


4- 1(T 5 


5.2 • 1(T 5 


2 • 1(T 5 



Table 2: Error values for the parameters given in table [T] of the fixed point solutions shown in fig. 



7.4. A„ s , A„ n and A a are the a-series error, the numerical error and the total error on the solution 
point in the a-plane. A& iS , Afe jn and Aj, have the analogous meaning for the solution point in the 
6-plane and Aoo )S , Aoo in and A ra for the asymptotic parameters A, B,C. All errors are relative 
errors. 

7.4 Structure of the a-plane 



From the configuration of lines in the a-plane on the right hand side in fig. 7.5 it might seem that 
the point where the brown and the olive line from J_ intersect the black solution line from I + on 
the singular line 73 represents a fixed point function f(R) defined on [R-, 00). As we will see now 



this is not the case. In performing the Taylor expansion (7.3) around zero we obtain expressions 



a n (ao, ai), n = 2, 3, ... such as (4.13), relating the higher coefficients to the first two. As is already 
the case for a2 these expressions have denominators that vanish along lines in the a-plane which are 



shown as dashed black lines in fig. 7.2 together with the solution lines. We have plotted the first 
four singular lines ji corresponding to singular denominators of the coefficients ai for % = 2, 3, 4, 5. 
For example, the singular curve 72 is given by 

3a (64vr 2 a + l) 



ai = 72(00) 



(7.14) 



192vr 2 ao - 7 

Each coefficient a n +i receives an additional factor in its denominator compared to a n and so we 
get a series of singular curves. The next one is 73, 

3a (256vr 2 a + 5) 



73 (ao) 



768vr 2 a - 25 
50 



(7.15) 



and similar expressions hold for all other singular curves. 

It is indeed straightforward to derive a general expression for the n-th singular curve j n . Sub- 



stituting the Taylor series (7.3) into the normal form of (3.2) we get 7 n by setting the coefficient 



of a n at order R n 2 to zero, n = 2, 3, The relevant terms can be picked out by hand from the 

normal form and one obtains the following condition: 



3 4 

-(ai + ao) + 7^ — ^ (I927r 2 (ai + a )a - 7a i + 3a ) = 



(7.16) 



(n-3)P ' (n-2)! 

This expression is also valid for n = 2 if we adopt the convention to keep only the second term in 
that case. If we solve this condition for a\ we obtain a general expression for j n . 



Let us now carefully look at 73. The constraint equation we get from plugging (7.3) into (3.2) 
at first order in R is given by 



9 ((— 7687r 2 ao + 25) a\ — 7687r 2 aQ — 15ao) 03 + 0220-2 + a2ia 2 + «20 = 0, 



where a<2 is known from zeroth order as in (4.13) and the a« depend on ao and a\ only: 



a 22 = 54 (-128vr 2 a + 3) 

a 2 i = 6 (-576vr 2 a 2 + 3 (-640vr 2 a + 3) a x - 640vr 2 a 2 1 - 7a ) 
a 20 = 8 (-144vr 2 a? - (4807r 2 a + l) a\ - 6 (48vr 2 a + l) oqOi) 



(7.17) 



(7.18) 
(7.19) 
(7.20) 



As explained before, the singular curve 73 is then simply given by setting the coefficient of 03 in 



(7.17) to zero. For any solution f(R) represented by a point on 73 we must also have 



«22a 2 + «2ia2 + «20 = 0. 



(7.21) 



Using (4.13), we can determine a discrete set of solutions for ao and a\ from the last equation 



together with the condition that we are on 73. Solving for them and calculating a 2 gives 



P3 



(a , ai) = (2.6589 • 10" 4 , 1.9684 • 10~ 4 ) , a 2 = 1.9913 • 10" 4 . 



(7.22) 



There are actually five possible solutions for the present case but we have singled out the unique 



non-zero solution corresponding to the point ps where the solution lines cross 73 in fig. 7.5 (right). 
All other solutions have ao > 7 • 10~ 3 and are thus far beyond the range we have developed the 
solution lines in. This entails that any solution line crossing 73 has to do so at discrete points 



such as p%. And indeed, this is exactly what we can see in fig. 7.5 all three solution lines, the 
two from I_ and the one from I + intersect 73 at p^, the only point on 73 in the plot which can 
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represent a solution. At this stage, the third coefficient a 3 is left undetermined at p%, and all higher 
coefficients 04, 05 . . . will be unique functions of 03. The parameter counting arguments now tell us 
that we expect a discrete set of solutions on 7 + and another discrete set of solutions on 7_ . This is 
a heuristic argument why we should not expect to find a solution at p% which is valid on 7_ U 7 + 
as these two discrete sets will in general not overlap. 

There is a way to determine the unique Taylor expansion at zero of any fixed point solution 
located on a solution line running through p% which proceeds without determining the free pa- 
rameter 03 by numerical integration. We should expect that the initial conditions at zero depend 



continuously on the initial conditions at We can then exploit (7.17) once more to determine 



03. The condition of continuity implies that if we solve (7.17) for 03 and take the limit as OQ,a\ 



approach P3 we will get a finite value 18 We can then proceed to obtain a unique Taylor expansion 
at P3 by exploiting the higher order constraints by which the coefficients 04,05, ... are fixed. Since 
we are taking limits of an expression depending on two variables in which both numerator and 
denominator become zero, it is not surprising to find the limit value depends on the direction along 
which we approach p^. If m denotes the slope of the straight line along which we take the limit, 
the function 

mi—)- lim 0,3(0,0, ai,m) (7.23) 

(ao,ai)-s>p3 

depends smoothly on m. Thus there is no reason to expect that the limits along the solution lines 
from 7 + and 7_ have to coincide. In fact, computing numerically the values of the slopes in each 
case and taking the limit, one finds the three different values 



A3H 



-1.297 • 10~ 4 , a|_ = 9.518 • 1(T 4 and a%_ = 7.239 • 1(T 4 (7.24) 



for the solution valid on 7+, the brown solution and the olive solution on 7_. 

What we have illustrated here in the case of the singular curve 73 holds analogously for all the 
other singular lines as well. This implies that the singular lines furnish a certain structure on the 
a-plane since solution lines are not free to cross them at any point. On each singular line there are 
only a few points that can be exactly determined where a solution line can go through. In general, 
the precise Taylor expansion at zero of the solution located at one of these points will depend on 
the slope with which the solution line crosses the singular line. 



Let us finally comment on the reason why it is imperative to use (7.9) with the behaviour shown 



in fig. 7.1 as a criterion to decide if we count a point as a solution point. Consider searching for 



3 whereas in the neighbourhood of P3 the same process for any point p € 73, will produce a divergent expression. 
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fixed point solutions on I + by integrating from R+ along the lines described in sec. 7.1, exploiting 



the corresponding form of (7.9) at zero instead of R+. Furthermore, let us first focus on the range 



ao > of the a-plane in fig. 7.9 , where we have plotted the first nine singular lines and the solution 
lines for J + . 

According to its definition the quantity Sf so \ can become arbitrarily small as we approach 
a solution line and shows singular behaviour as we cross a singular line of the a-series. More 
precisely, since the n-th singular line is only relevant for coefficients Qfc(ao, «i) with k >n and we 
are working with an a-series of order 5, 5f so \ will only show singular behaviour close to the lines 
72, . . . , 75. Away from these two special cases <5/ so i has a typical size that varies across the a-plane 



(as long as we stay above 72 in fig. 7.9 numerical integration from R + to zero is generically possible 
despite the threat posed by moveable singularities). Speaking in orders of magnitude, 6f so \ starts 
at 10 -3 between 72 and 73 and drops to 10~ 8 in the area between 74 and 75. If we now take into 
account the next order in the Taylor expansion around zero, i.e. we include the term with coefficient 
ag(ao,ai), we find that its size is always several orders of magnitude smaller than 5f so \ in these 
regions and thus will have no effect on existing solution lines. This changes once we consider the 
region above 75 towards 76 and beyond. Here, the critical quantity Sf so \ becomes comparable in 
size to the 6th order correction of the a-series, both evaluating to around 10 -14 as we approach 
76- When searching for solutions in a region above the singular curve j n it is therefore essential to 
use a Taylor series at zero of order at least n + 1. For example, using a Taylor series of order 6 at 
zero we found a first solution point of a solution line situated between 75 and 76 in the region of 
positive ao and we would expect that more solution lines can be found continuing this pattern. 

Coming back to why we should not just require |£/ so i| to be smaller than some upper bound 
for a point to be regarded a solution point, it is clear from this discussion that there is no sensible 
way to find such a bound. Moreover, the value of |5/ so i| depends sensitively on eq. For example, 



with a 5th order Taylor series at zero and £0 = 5 • 10~ 5 a scan of the type in fig. 7.1 between 75 
and 76 at ao = 7 • 10 -4 produces values of <5/ so i between 10~ 12 and 10~ 14 but Sf so \ never changes 
sign. Decreasing £0 in several steps to £0 = 10 -6 , the same scan shows correspondingly decreasing 
values of |<5/ so i| reaching 1CP 18 to 10 -21 in the last step. During none of these steps there is a sign 
change of 5f so \ and thus no evidence for a solution line. This again shows that just imposing an 
upper bound on |5/ so i| cannot be sufficient. It also confirms the necessity of the 6th order term in 
the a-series if one intends to find solutions between these two singular lines. 



What we have described here in the region ao > of the a-plane in fig. 7.9 holds analogously in 
the region of negative ao- Indeed, the bottom part of the large solution line on this side has been 
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Figure 7.9: The solution lines (black) of fig. 7.2 in the a-plane together with the singular lines (red). 



From bottom to top in the range of negative ao the singular lines are 72, 73, ... , 710 stemming from 



the coefficients 02, as, ... , aio of the a-series (7.3). 



developed with a Taylor series of order 6 at zero from the beginning. 

One might also wonder if there are similar singular lines in the 6-plane arising from the expansion 



(7.2). An investigation into the expressions for its coefficients indeed shows that such singular lines 
are present. However, the solution lines we have found in the 6-plane are all a safe distance away 
from them and we therefore do not expect a similar picture. 



8 Discussion and Conclusions 

In this section we highlight the main technical developments, then the main discoveries and their 
physical interpretation, and finally discuss the extent to which we should trust these discoveries 
given the methods used to derive the f(R) approximations in the first place. 

Our main overall conclusion is a technical one: the parameter-counting arguments, introduced in 
scalar field theory at the LPA and higher levels of derivative expansion |19| - |23| and adapted here to 
this context, remain a very powerful way to determine what to expect from the equations in terms of 
fixed points and perturbations. The counting argument is the statement that any restrictions to the 
dimension of the parameter space arising from fixed singular points and the asymptotic behaviour, 
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act independently, and therefore, if there are any global (or partial) solutions, the dimensionality of 
their parameter space can be determined by counting up the number of such restrictions. Although 
the application in this paper is to the asymptotic safety programme in quantum gravity, this simple 
counting argument should be applicable far more generally, providing crucial insight whenever one 
is searching for smooth solutions to highly non-linear renormalisation-group-type equations with 



little or no inherent symmetry (c/. sec. 4.2). 

That this insight is crucial is underlined by our numerical investigation in sec. [7| The basic 
method of attack is as follows: series expansions are developed analytically around the fixed singular 
points and around infinite field, and used to 'patch' onto numerical solutions in between. Apparent 
violations of the counting rules arise quite often during numerical investigations of the resulting 
parameter space of partial solutions. If one trusts the rules, this motivates further careful numerical 



and analytical studies at these points. As explained in sec. 7.4 we found a whole region in which 
S /soi, as given by the corresponding version of (7.9) at zero, fell to less than 10 -14 . By decreasing 
£o this fell to 10 -21 , but never changed sign. We insist that a sign change is required in order to 
identify a solution. This phenomenon led to the recognition of the importance of higher terms in 
the a-series and the corresponding curves in the a-plane where these coefficients diverge (except at 
special points). As a result, we uncovered a beautiful structure which is effectively responsible for 
the lines of solutions and guides them through the a-plane, cf. fig. |7.9| 

Important for the application of the counting argument is a thorough understanding of the 
analytic asymptotic expansion of the fixed point solutions and their eigen-perturbations. We showed 
how to check that the full parameter space of asymptotic behaviour has been uncovered, and how to 
develop systematically the higher orders, even in cases such as the current one where the differential 
equations and the asymptotic series are highly involved. Again, these technical developments are 
valid much more generally, thus we expect them to be important in applications other than to 
asymptotically safe quantum gravity. 

Now we turn to the discoveries we made by applying these methods to the f(R) approximations 
for asymptotically safe quantum gravity as introduced in refs. |15|4l7|. Using the counting argument 



we showed that the f(R) approximations derived in refs. 15 , 16] have no global smooth solutions, 
and thus are unsatisfactory. 



The reason can be traced to fixed singular points introduced via their choice of cutoff 17 



In ref. 17 , the cutoff was chosen so that these were avoided, however three fixed singular points 
remain: R = and R = R±. We showed that the asymptotic series for large R developed in 
ref. 117] was only part of the full parameter space of asymptotic solutions. The full asymptotic 
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series solution provides no restriction on the dimension of the parameter space but only restricts 



the three parameters, A,B,C, to lie inside a cone (5.17). 

The f(R) approximations are derived on four-spheres, and thus we have < R < oo. Actually 
the point R = must also be included to satisfy the quasi-locality requirement of the Wilsonian 
renormalisation group as discussed in sec. [3j Since only the R = and R = R + singular points 
are operative, the counting arguments show that continuous lines of fixed points should exist. We 
confirmed this in sec. [7j by uncovering five separate lines of fixed points, and evidence for a sixth 
line. As we saw, there are very possibly more lines, but uncovering these will require further 
increases in the number of significant digits kept during the numerical calculations. 

A fixed point solution on one of these lines is determined for example by providing the initial 
conditions /(0) = ao and /'(0) = a\. We showed how the existence and position of these lines in 
(ao,oi) space can be understood from the singularity structure of the Taylor series coefficients at 
R = 0. We have uncovered lines with a\ < as well as lines with at > 0, corresponding to both 
possible signs of the coefficient for the y/gR term in Minkowski signature. Similarly, the cosmological 
constant term, corresponding to ao, is negative on some of the lines of fixed points we analysed and 
positive on others. However, the sign of —a\ and ao have no direct connection to the sign of the 
physical Newton's gravitational constant G and cosmological constant A that we would predict from 
this theory. To compute these latter quantities would require adding all the (marginally) relevant 
perturbations about the fixed point and integrating out along the corresponding renormalised 
trajectory, taking the limit as k — > in such a way as to obtain an effective f(R) parametrised 
by the finite (marginally) relevant couplings. The predicted Newton's constant and cosmological 
constant then follow from its Taylor expansion (see for example the corresponding situation for 
scalar field theory in ref. |22|). 

Of more importance is the fact that we have found lines of solutions with either sign for A. 
Only the solutions with A > have a chance to be physically sensible since these correspond to 
fixed point solutions f(R) that are bounded below (in both Euclidean and Minkowski signature) 
and which in unsealed physical variables having completed the functional integration (i.e. in the 
k — > limit) result in a positive f(R) = AR? . (Note that this does not mean that the fixed points 
correspond to an R 2 type theory: infinitely many equally important interactions have been sent to 
zero in this limit.) 

We developed the full asymptotic expansion for infinitesimal perturbations around these fixed 
points in sec. |6.1|and furthermore analysed the asymptotic behaviour of finite perturbations in sec. 



6.3 Unlike in the LPA for scalar field theory, we saw no reason to exclude some of this asymptotic 
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behaviour. This is not surprising, given that here the asymptotic expansion of the fixed point itself 
does not restrict the dimension of the parameter space. In this case the counting arguments show 
that each fixed point supports a continuous infinity of operators, and operator dimensions are not 
quantised. Therefore we are forced to conclude that all these fixed points are unsatisfactory. 

We saw that the situation can be resolved if we also include consistently spaces with —oo<R< 
0, as we are required to do, since quantum gravity should also make physical sense for such spaces. 
We therefore searched for smooth global solutions over the full range — oo < R < oo. The counting 
arguments now tell us that any such fixed points are discrete in number and support a quantised 
eigen-operator spectrum. However, after developing a detailed understanding of the structure of the 
parameter space, and performing an extensive numerical search, we failed to uncover any solutions 
in this case. The closest we came to finding such solutions numerically, were partial solutions which 
do not extend to R = — oo but which are defined on the interval [R~, oo). 

To what extent can we trust all these discoveries as uncovering the real physics of the situation 
rather than artefacts of the approximations and methods used? Although, as we pointed out in 
the introduction, truncation to general functions f(R) is still severe, breaking in an uncontrolled 
way both scheme invariance and the underlying diffeomorphism invariance of the quantum field 
(expressed through modified BRS Ward identities [2|), and further followed by smoothing approx- 
imations, we will see that some of the effects clearly have their origin in some generic qualitative 
features and whose presence can be justified on physical grounds. 

As we have sought to demonstrate, our methods for analysing the resulting differential equations 
are powerful and exhaustive and we feel confident therefore that the parameter space has been fully 
explored, except in the following respects. Firstly, as already mentioned, there are almost surely 
further lines of fixed points valid in U loo whose development requires additional terms in the 
a-series. Secondly, at the expense of time-consuming numerical integration, it would be possible to 
extend some of the known solution lines further into the region of larger values of ao where their 
exact behaviour is still to be determined. 

The fact that we need to consider consistently also R < surely reflects the real physics of the 
situation. However, we are now using the equations outside the four-sphere region in which they 
were derived. As we noted in sec. [TJ really we should use an f(R) approximation that consistently 
takes into account throughout its derivation, spaces in which R can be negative. This task is beyond 
the scope of the present paper. 

Even if we restrict ourselves to R > 0, we uncovered many novel features and these are conse- 
quences of the fixed singular points R = 0, i?+ and the novel asymptotic behaviour. 
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As discussed at the end of sec. [2l the R = singular point arises as a consequence of choosing 



cutoffs that amount to modifying the Laplacian, cf. (2.6), thus forcing the cutoff to depend on 
/". This also turns the fixed point equation from a second order ordinary differential equation, 



as would have been expected from the structure of (2.1), into a third order one. The parameter 
counting is thus unaltered, suggesting that the R = singular point is not an unwanted artefact 
but is in fact necessary, reflecting the true physics of the situation [l~7| . 

In the LPA we find that the leading asymptotic behaviour is determined by dimensions, and 



this is because quantum fluctuations can be neglected at large field ip 1 19 , 23 . Here we find as in 



ref. 17 . that the large field behaviour of the fixed point solution is again as expected by dimensions. 



But the reason is more subtle, as discussed below (5.12). This is reflected in the fact that for the 



perturbations, the asymptotic behaviour is not determined by dimensions - see below (6.6). 

Recalling that we multiplied through by R 2 in (2.12), we see that the quantum fluctuations for 
the fixed point solution are 0(1) in this regime, whilst for perturbations they vanish as 0(1/ R). In 
scalar field theory, the quantum fluctuations in the sharp cutoff case are also 0(1) (up to logarithmic 
corrections) fl9| whilst for smooth cutoffs they vanish as ip — > oo, see e.g. [20] . The crucial difference 



is that here the left hand side of (2.12) is multiplied by the volume which is vanishing in the large 
field limit, as 1/R 2 , making the classical contribution subdominant in the large field regime. This 
points to a fundamental physical reason for the novel asymptotic behaviour: at the same time as 
R — > oo, where one might naively expect to be able to neglect quantum corrections in the Wilsonian 
renormalisation group, the volume of the space is vanishing and this forces quantum fluctuations 
to remain important. It would be desirable to strengthen this conclusion through a more thorough 
investigation of the asymptotic dependence of quantum fluctuations on choice of cutoff. 

This is all the more important since the physical scalar quantum fluctuations determine even 
the leading asymptotic behaviour of the perturbations, as we saw from the discussion surrounding 



(6.6). As we noted below (5.12), the leading R 2 behaviour for the fixed point itself comes about 
by balancing the physical scalar quantum corrections against the non-physical scalar quantum 
corrections. It is tempting to speculate that these effects are somehow a reflection of the instability 
present in the conformal mode, which at first sight leaves no imprint on the exact renormalisation 
group [2j. 

The final key feature was the presence of the fixed singular point at R = As discussed 
in ref. |17| and reviewed at the end of sec. § this arises from the failure to integrate out the 
constant physical scalar mode. More study seems called for in order to establish whether this 
is really a generic feature of cutoff choices that are otherwise acceptable [17], and to understand 
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the consequences, both from the renormalisation group point view and more generally, of leaving 
behind such an unintegrated mode. 

Note that intuitively it is clear that the underlying diffeomorphism invariance (as expressed 
through modified Ward identities) does not allow complete freedom in choosing separate cutoffs for 



the component fields in (2.3) and the ghosts and auxiliary fields, otherwise we would fail completely 
to recover the fact that the components must all effectively combine back into the quantum field 
hau- However the extent to which the cutoff choices are constrained, and the consequences of not 
obeying such constraints, has yet to be explored in the literature. 



Returning to the earlier forms of the f(R) approximation [15,16 , we have already seen that 
their choice of cutoff resulted in fixed singular points which in turn ruled out the existence of 
any global smooth fixed point solutions. Similar, but /(-R)-dependent, single poles are responsible 
for the unphysical termination of certain renormalised trajectories in polynomial truncations (e.g. 
termination of type Ilia trajectories as described in ref. [l0l[T5]). It has been conjectured that 



these terminations are an artefact of the truncations used 15, 16 and indeed it has been shown 



that they can be avoided in a suitably chosen larger space of operators, but this is at the expense 



of introducing non-polynomial functions of R which violate quasi-locality 15 . In ref. 16 , it was 



noted that the termination of these trajectories can be avoided if so-called type III cutoffs are used. 



This mechanism bears similarities to the avoidance of fixed singularities in (2.12) by use of a cutoff 



of type II (see the discussion below (2.11) and ref. [17]). 

Clearly this all points to the fact that care needs to be exercised over the choice of cutoff. 
Apparently sensible choices of cutoff can result in a break-down of the Wilsonian renormalisation 
group, perhaps only in certain regions of parameter space, and this is evidenced by unphysical 
behaviour. In fact such a situation is well established in simpler systems and in some cases rigorously 
proved [36] . 

We see from these investigations that unsurprisingly the form of the cutoff plays a crucial role 
in the nature of the solutions. In any case, in view of the novelty and number of unexpected 
features that we have uncovered, more investigations using different cutoffs and different types of 
approximation at this LPA-like level of truncation, or beyond, are evidently called for. 
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Appendix. An early version of the fixed point equation of ref. [ 17| 



At an earlier stage of this work, before we became aware of an error in the first version of ref. 17 



we had carried out an analysis of the uncorrected version of the fixed point equation (3.2). This 
first form of the equation is obtained by deleting the term 36(i? 2 + 2)(/' + 6/") in the contribution 



of the physical scalar modes in (2.16) and the corresponding term in (p. 4b. This change does not 



affect the technical advances in this work such as the validity of the parameter counting method, 
the structure of the asymptotic expansion and eigen-perturbations, or the conclusions drawn from 
these. The main features of the set of solutions of the fixed point equation we had uncovered also 
resemble the results presented in sec. [7j We describe these solutions here to illustrate the general 



applicability of the methods used in this work for analysing such equations as (3.2). 



In total we found four solution lines in the range i+ U 1^ and one solution line on the interval 



J_. 
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Figure A.l: The four lines of solutions in the 6-plane. Different solution lines are represented 
consistently in different colours (from left to right these are black, red, green and blue). The bold 
part of each line distinguishes solutions valid for all R > from those valid only on /_)_. 



They are shown in fig. A.l and A. 2 One main difference to the results in sec. [7] is the existence 
of a compact subset for each solution line for I + containing only those points that define fixed 
point solutions valid on I + U 1^ with asymptotic parameters A, B, C satisfying a modified version 
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Figure A. 2: The four lines of solutions in the a-plane corresponding to the Mines in fig. A.l (from 



bottom to top in black, red, green and blue) together with the only solution line for J_ in brown. 
The solutions valid for all R > are again given by the bold parts of the lines. Along each of the 
black dashed lines a particular coefficient in the Taylor series around zero becomes singular. 



of (5.17), 



^A 2 >B 2 + C 2 . 
20 



(A.2) 



With regard to solutions for negative values of R the only solution line we were able to find is also 



shown in fig. A.2 The corresponding solution line in the /3-plane is displayed in fig. A. 3 It does 
not contain a bold part which means that no solutions valid for all R < were found. Like the 
solution lines in the /3- and 6-plane of sec. [7j the solution line in the /3-plane is mostly linear with 
only small non-linear variations superimposed. In this analysis we did not use both roots of the 
quadratic constraint for the second coefficient of the b- and /3-series which may be why the solution 
lines end at the origin of their respective plane. 

Just as in sec. [7j we also have singular lines associated with the Taylor expansion around zero 



shown as black dashed lines in fig. A.2 Furthermore, all the conclusions drawn in connection with 



singular lines in sec. 7.4 hold in the present situation. In particular, solution lines are allowed to 
cross singular lines at only a finite number of points on them as is beautifully illustrated by the 



brown solution line for /_ in fig. A.2 cutting through the singular lines precisely where the solution 



lines for I + do. Similarly, we do not obtain fixed point solutions valid on I_ U i+ at these points 
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o o.oi 0.02 0.03 o.'o4 ~o am ao2 ao3 ' ao7 

Figure A. 3: The (only) solution line in the /3-plane is shown on the left together with its nonlinear 
variations on the right. The straight line m ■ + c that has been subtracted from the line on the 
left has slope m = -0.2312 and c = 1.14 • 1CT 5 . 

since the value of the coefficient in the a-series associated with the singular line depends on the 
slope with which the solution lines cross the singular line. 
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